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ABSTRACT 

We present evidence that the incidence of active galactic nuclei (AGNs) and the distribution of their 
accretion rates do not depend on the stellar masses of their host galaxies, contrary to previous studies. 
We use hard (2-10 keV) X-ray data from three extragalactic fields (XMM-LSS, COSMOS and ELAIS- 
Sl) with redshifts from the Prism Multi-object Survey to identify 242 AGNs with L 2 _ w koV = 10 42 " 44 
erg s~ x within a parent sample of ^25,000 galaxies at 0.2 < z < 1.0 over ~ 3.4 deg 2 and to i ~ 23. 
We find that although the fraction of galaxies hosting an AGN at fixed X-ray luminosity rises strongly 
with stellar mass, the distribution of X-ray luminosities is independent of mass. Furthermore, we 
show that the probability that a galaxy will host an AGN can be defined by a universal Eddington 
ratio distribution that is independent of the host galaxy stellar mass and has a power-law shape with 
slope —0.65. These results demonstrate that AGNs are prevalent at all stellar masses in the range 
9.5 < XogM.*/ M.q < 12 and that the same physical processes regulate AGN activity in all galaxies 
in this stellar mass range. While a higher AGN fraction may be observed in massive galaxies, this is 
a selection effect related to the underlying Eddington ratio distribution. We also find that the AGN 
fraction drops rapidly between z ~ 1 and the present day and is moderately enhanced (factor^ 2) 
in galaxies with blue or green optical colors. Consequently, while AGN activity and star formation 
appear to be globally correlated, we do not find evidence that the presence of an AGN is related to 
the quenching of star formation or the color transformation of galaxies. 
Subject headings: galaxies: active - galaxies: evolution - X-rays: galaxies 



1. INTRODUCTION 

Supermassive black holes (SMBHs) are now thought 
to reside at the centers of most, if not all, galaxies 
with central stellar bulges (IKormendv fe Richstondll995l : 
iKormendv fe Gebhardtil2001[ ). These galaxies must have 
undergone periods of intense accretion activity to build 
their black holes, during which they would be observed 
as active galactic nuclei (AGNs). The mass of the 
SMBH is kno wn to be tightly corr elated with both 
the lumi nosity dMagorrian et al.lll99cf) and velocity dis- 
persion dFerrarese fc Merrittil2000t IGebhardt et al.ll2000t 
iTremaine et al.l I2002D of the stellar bulge, indicating 
that the growth of the SMBH and the processes that 
build up the stellar mass of the host galaxy must be 
related, although it is unclear whether a ca usal connec- 
tion between these pro cesses is required (see lPeng]|2007l : 
Uahnke fc M accio 2011). In addition, the star formation 
rate density as a function of redshift is found to correlate 
with the black hole accretion rate density out to high 
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redshifts fe.g.. iSilverman et al.ll2008l: lAird et~aTl 120101) , 
indicating that these processes may be controlled by a 
common mechanism. Indeed, a number of studies have 
proposed that AGNs provide feedback that quenches star 
formation in their host galaxies and is crucial to con- 
trol the transformation of blue, star- forming galaxies to 
red, passively evolvin g systems (e.g..lCroton et al.l 120061 : 
Hopkins et all l2007at ISchawinski et al.ll2007t ). although 
there remains a lack of direct observational evidence of 
this process. Nevertheless, what triggers and regulates 
AGN activity and how this relates to the co-evolution of 
AGNs and their host galaxies remain key, open questions 
in extragalactic astrophysics. 

To understand how AGNs and galaxies co-evolve we 
must study the types of galaxies that tend to host AGNs. 
Galaxies can be separated into a variety of classifica- 
tions, many of which have significant overlap. Many 
studies have shown that the positions of galaxies in the 
optical color-magnitude diagram are strongly bimodal, 
falling into two distinct regions associated with passively 
evolving galaxies on the "red sequen ce" and star-forming 
galaxies in the "blue cloud" (e.g. , iBlanton et al.l 120031 : 
iBaldrv et all2004t iBell et al.ll2004f ): a narrow "green val- 
ley" separates these two populations. The total stel- 
lar mass of a galaxy, Aa,, is also a key parameter that 
tracks the net result of star formation activity through- 
out a galaxy's lifetime, as well as the hierarchical build 
up of galaxies through mergers. While both the red 
and blue galaxy populations span a wide range of stel- 
lar masses, red galaxies tend to dominate the popula- 
tion at higher stellar masses, whereas blue galaxies gen- 
erally have lower stellar masses. Over cosmic time the 
total combined mass of galaxies on the red sequence is 
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found to increase, indicating that star formation must 
be quenched in blue galaxies, which subsequently allows 
their stellar populations to evolve pas sively and trans- 
fers the galaxies to the red sequen ce (jBell et al.1 120041 : 
iFaber et al.ll2007HBrown et al.ll2007| ). Morphology is also 
closely correlated with galaxy color and stellar mass — 
massive, red galaxies are generally early-type, bulge- 
dominated systems, whereas blue galaxies tend to have 
late-t ype morphologies with subst antial disk components 
(e.g.. iBlanton fc Moustakasll2009ft . 

Which of these galaxies host AGNs? Given the corre- 
lation between the ov erall AGN accretion and star for- 
matio n rate densities ()Silverman et al.l 120081 : lAird et all 
12010ft . we might expect to find AGNs predominantly 
within blue star-forming galaxies, which tend to have 
low stellar masses. A number of studies show that 
this is not the case — the fraction of galaxies hosting 
AGNs increases in higher stellar mass, bulge-dominated 
galaxi es — although the s ubject remains under much de- 
bate. Kauffm ann et al.l (|2003al ) found that narrow-line 
AGNs (identified by their emission-line ratios) are pref- 
erentially in massive, bulge-dominated galaxies (A-f*> 
1Q W M Q ) in the local universe. iKauffmann et al.l (|2003aft 
also found evidence for recent star formation in the 
host galaxies of the more luminous AGNs, although 
the host galaxies of the lower luminosity AGNs gener- 
ally had older stellar populations, similar to mo st early- 
type galaxies (see also Kauffm ann et al.l l2007f ). The 
fraction of galaxies that host radio-loud AGNs is even 
more strongly dependent on stellar mass, rising rapidly 
from < 0.01 % at M^ 10j° M Q to > 30% at M*> 
5 x 10 n A4 g (lBest et al.l l2005h : it has long been estab- 
lished that radio-loud AGNs a re primarily found in 
massive, elliptical g alaxies (e.g.. iMatthews et al.l [1964J : 
iDunlopet aLll2003ft . iSchawinski et al.l (j2007ft studied the 
optical colors of elliptical galaxies in the local universe 
and found that those hosting narrow-line AGNs were 
generally on the red sequence; these results suggest that 
star formation and AGN activity do not occur at the 
same time, which could be explained by the onset of 
AGN activity (possibly during an early, highly obscured 
phase of the AGN lifetime) being responsible for quench- 
ing star formation activity within galaxies and transfer- 
ring them from the b lue cloud to the red sequence (e.g., 
iHopkins et alJl2006aI ). However, direct observational 
evidence of the AGN feedback process is lacking. A de- 
tailed study of the prevalence of AGNs as a function of 
stellar mass and galaxy color is required to shed light 
on the relationship between the growth of galaxies, the 
progress of star formation, and the accretion activity of 
the central SMBH. 

Deep X-ray surveys have enabled efficient identifica- 
tion of obscured and low-luminosity AGNs where the 
host galaxy dominates the optical light, allowing stud- 
ies of the host galaxy prope rties to be push ed out to 
much higher redshifts (z < 4; iXue et all 120101 ) . Various 
studies have shown that the fraction of galaxies hosting 
X-ray AGNs rises strongly wi th increasing M.* at all red- 
shifts from z ~ 0.4 - 4 (e.g.. lAlonsq-Herrero et"aLll2008l: 
iBrusaet al.ll2009bl: IXue et all 12010ft. similar to the op- 
tical narrow-line AGNs id entified in the local universe 
(jKauffmann et al.l l2003ah . iNandra et all (|2007l ) investi- 
gated the positions of X-ray AGN hosts in the optical 



color-magnitude diagram at z ~ 1 and found that X-ray- 
selected AGNs are predominantly in luminous galaxies 
but avoid the strongest star-forming galaxies and thus 
lie in the reddest part of the blue cloud, in the green 
valley, or on the red sequence. These results indicate 
that AGN activity may be associated with the transfor- 
mation of blue galaxies and the quenching of star for- 
mation, indicating that AGN feedback could potentially 
pl ay a key rol e in th is process. Similar results were found 
bv lCoil et ail (|2009l ) for a larger sample of X-ray AGNs, 
and t he same pattern appears to hold at lower red - 
shifts dHickox et al.ll2009t iGeorgakakis fc Nandr a! 12011ft 
However. ISilverman et all (|2009aft and IXue et al.[ (j2010ft 
have demonstrated the importance of st ellar-mass selec- 
tion eff ects, which may bias prior work. ISilverman et al.l 
(2009a) found that X-ray AGN activity was associated 
with ongoing star formatio n, and thus the A GN fraction 
is higher in the blue cloud. IXue et all (j2010ft found that 
the AGN fraction rises strongly with increasing A4*, as 
described above, but when considering samples matched 
in stellar mass there was no evidence for specific cluster- 
ing of AGN host g alaxies on the color-magnitude dia- 
gram. Cardamone et al.1 (|2010ft also studied the relation 
of AGN host galaxy masses and colors; they found that 
the colors of AGN host galaxies of a given stellar mass 
followed the same bimodality as the non-AGN galaxy 
population, after correcting for the effects of dust. Thus, 
despite much recent progress, it remains unclear as to 
whether there is an association between AGN activity 
and the color transformation of galaxies, and how this 
relates to the predominance of AGNs in massive galax- 
ies. 

Further insight into the processes that trigger AGNs 
and how these relate to the evolution of galaxies may 
be possible by directly probing the distribution of AGN 
activity, the active growth of SMBHs via accretion that 
produces a luminous AGN. A number of studies have 
constrained the shape and evolution of the AGN lumi- 
nosity function, wh ich traces AGN activity, over a wide 
redshift range (e.g. , lUeda et al.l [2003: Eb rero et al.ll2009b 
lYencho et all 1200ft lAird et alJ 1201(1 : lAssef et al.N20Tlft . 
The luminosity function has a steep double power-law 
shape with a break at a characteristic luminosity, L». 
The rapid decline in the AGN accretion rate density since 
z ~ 1 is a result of a downward shift in L» for the over- 
all population, but the l uminosity function retains es- 
sentially the same shap e (|Barger et al.1120 05; Air d et aLl 
120101 : lAssef et al.l 12011ft . This indicates that the pro- 
cesses responsible for triggering AGNs, fueling black hole 
growth and controlling the lifetimes and duty cycles of 
AGNs are essentially the same but are moving towards 
generally less luminous systems between z ~ 1 and the 
present day. However, few observational studies have at- 
tempted to simultaneously study the distribution of AGN 
activity as traced by the AGN luminosity and the preva- 
lence of AGNs as a function of the host galaxy properties 
such as color and stellar mass. 

Understanding both the predominance of AGNs in 
massive galaxies and any correlation with host galaxy 
colors is vital to reveal the role and effect of AGNs on 
the evolution of galaxies. In addition, to constrain the 
processes that trigger AGNs and drive black hole growth 
we must track the evolution of AGN accretion rates and 
the masses of their central SMBHs, although direct mea- 



PRIMUS: The dependence of AGN accretion on host stellar mass and color 



surement of these quantities is extremely difficult, if not 
impossible for large samples of distant AGNs and galax- 
ies. In this paper, we use the X-ray luminosity as a tracer 
of the AGN accretion rate and investigate the relation- 
ship with the total stellar mass and colors of galaxies. 
While the relationship between these observed proper- 
ties and the underlying physical parameters remains un- 
certain, careful investigation including consideration of 
observational biases and selection effects is vital to re- 
veal and constrain the connection between the growth of 
galaxies and their SMBHs. Our study focuses on inter- 
mediate redshifts {z < 1), a key time period of the uni- 
verse over which star formation and AGN activity have 
rapidly declined from the peak at z ~ 1 — 2. We use the 
PRIsm MUlti-object Survey (PRIMUS: ICoil et all 12011 
Cool et al, in preparation), a recently completed faint- 
galaxy survey that used a prism to obtain low-resolution 
spectra and determine redshifts for ~ 120,000 galaxies 
over ~ 9 deg 2 with z ~ 0.2 — 1.2, providing the largest 
currently available spectroscopic sample of galaxies at 
these redshifts. X-ray data are available for > 3 deg 2 of 
PRIMUS from a variety of Chandra and XMM-Newton 
survey programs. We use these X-ray data to identify 
AGNs within the PRIMUS sample, enabling us to effi- 
ciently select AGNs down to low luminosities and with 
moderate obscuration. For such sources the optical light 
probes the host galaxies, and we can compare properties 
with the extremely large non-AGN galaxy sample pro- 
vided by PRIMUS. Section [5] briefly describes our data, 
providing details on the compilation of X-ray source lists 
and matching to PRIMUS sources. Section [3] outlines 
our derivation of stellar masses and Section @] describes 
in detail the selection of our AGN and galaxy samples. 
We investigate how the probability of a galaxy hosting an 
AGN is related to the stellar mass of the galaxy and the 
AGN X-ray luminosity in Section [5] for two wide redshift 
bins. In Section [5] we determine the redshift dependence 
of the AGN fraction and perform a more robust analy- 
sis of the stellar mass and luminosity dependence that 
accounts for the evolution with redshift via a maximum- 
likelihood fitting approach. In Section [7] we show that 
the probability that a galaxy hosts an AGN is primar- 
ily determined by a power-law distribution of Eddington 
ratios with a slope of —0.65 and a normalization that 
depends on redshift but does not depend on the stellar 
mass of the host galaxies. We investigate the depen- 
dence of the AGN fraction on galaxy color in Section [5] 
and find a weak enhancement of AGNs in galaxies with 
bluer colors. We discuss the implications of our results in 
Section [S] Section [TU] summarizes our results and overall 
conclusions. 

Throughout this paper we adopt a flat cosmology with 
51a = 0.7 and h — 0.7. All magnitudes are on the AB 
system. 

2. DATA 

2.1. PRIMUS 

PRIMUS is the largest faint-galaxy, intermediate- 
redshift survey performed to date and covers ~ 9 deg 2 
in a total of seven different well-studied multiwavelength 
fields. We used the IMACS instrument on the Magellan I 
Baade 6.5m telescope with a slitmask and low-dispersion 
prism to obtain low-resolution (A/AA ~ 40) spectra for 



~ 2000 objects with a single slitmask. We generally ob- 
served two slitmasks for each pointing, allowing us to 
assign slits to a statistically complete sample down to 
i « 23. At our faintest magnitudes (i ~ 22.5 — 23) we 
sparse-sample galaxies with a weight of 0.3 if an object 
has potential slit collisions; this sparse-sampling ensures 
that PRIMUS is not dominated in number by the faintest 
galaxies. We also applied a density-dependent sampling 
scheme to all objects with more than two nearby neigh- 
bors on the plane of the sky to reduce the effects of slit 
collisions in a well-controlled manner. Overall we were 
able to assign slits to ~ 80% of galaxies at i < 22.5. 
Unresolved point sources were only retained if their op- 
tical and/or IR colors (or an X-ray detection for the 
Subaru/ XMM-Newton Deep Survye, SXDS, area of the 
XMM-LSS field) indicated they may contain an AGN; 
however, no additional prioritization was given to such 
sources in the targeting or slitmask design. For full de- 
tail s of the survey de sign, targeting and a data summary 
see ICoil et all (12011 . 

To obtain redshifts we fit the low-resolution spectra 
and photometry with galaxy as well as BLAGN and stel- 
lar spectral templates. We supplement our low-resolution 
spectra with optical ground-based photometry, which 
is used in the redshift fitting as well as to derive K- 
corrections and stellar mass estimates (see Section [3] be- 
low). We classify the objects as stars, BLAGNs or galax- 
ies based on the \ 2 of the best fits. Our redshift success 
rate is ~ 90% at bright magnitudes (i < 21) and drops to 
~ 70% at i ~ 22.5, but has little dependence on galaxy 
color. Our final redshift catalog contains ~ 120,000 ro- 
bust redshifts with a precision of a z /(l+z) =0.005. (nor- 
malized median absolute deviation) For details of the 
data reduction, redshift fitting algorithm, redshift confi- 
dence and precision, and survey completeness see R. J. 
Cool et al. (in preparation). 

2.2. X-ray data 

X-ray data of varying depths from Chandra and/or 
XMM-Newton are available for five of our PRIMUS 
fields. For this paper, we have compiled publicly avail- 
able X-ray point source catalogs in three of these fields: 
XMM-LSS, COSMOS and ELAIS-S1. Further details of 
the catalogs in the individual fields are given in Sections 
|2~2~T1 12X21 and [2X3] below. We note that these X-ray 
catalogs were compiled after the PRIMUS observations 
had been completed. PRIMUS is a purely flux limited 
survey and the X-ray information had no bearing of the 
targeting of sources. 

In all our fields we h ave used the likelihood ratio 
matching t echniqu e (e.g., [S utherland & Saunder s] Il992t 
iCiliegi et a l. 2003; Brusa et al.l2007HLaird et al.ll2009l Uo 
robustly identify the counterparts of each X-ray source in 
the PRIMUS targeting opti cal photometry ( usually the 
i or R band, see Table 3 in ICoil et all 120111 ). The like- 
lihood ratio technique allows us to assign secure optical 
counterparts to the X-ray sources, accounting for the op- 
tical and X-ray positional uncertainties, the probability 
of a counterpart with a given magnitude, and the prob- 
ability of a spurious match. Candidate counterparts are 
identified within 5" of an X-ray source position. While 
a 5" search radius is significantly larger than the typical 
X-ray positional uncertainties for on-axis Chandra data, 
it is appropriate for XMM-Newton or off-axis Chandra 
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data. As the likelihood ratio takes account of the posi- 
tional uncertainties for individual sources, we can adopt 
a single search radius for all our fields and data. Where 
multiple counterparts exist we choose the match with the 
highest likelihood ratio; we further restrict our matches 
to "secure" counterparts with likelihood ratios > 0.5. 

We match these optical counterparts to sources actu- 
ally targeted by PRIMUS and inspected the PRIMUS 
spectra of every X-ray source by eye. These manual 
checks allow us to confirm the PRIMUS classifications 
and redshifts, identify spectra with clear features (such as 
broad lines) that may not be well-fitted by the automated 
procedure, select between multiple possible redshift so- 
lutions and identify catastrophic data issues. Based on 
these manual checks we altered the redshifts and clas- 
sifications of around 20% of our X-ray sources, mainly 
reclassifying sources as BLAGNs. 

After completing our manual checks, we compiled prior 
estimates of the redshifts of our sources based on high- 
resolution spectroscopy from other surveys. I n the COS- 
MOS field we use the zCOSMOS catalog ([Lilly et al.l 
:2009f) as well as redshifts from various spectroscopic 
campaigns following up XMM-COSMOS sources, as pre - 
sented in the recent compilation by Bru sa et al.l (2010). 
In the XMM-LSS we use the compilation of spectroscopic 
redshifts in t he UPS aree0, which i nclu des re dshifts from 
iSmail et all (12008ft . ISimpson et al.l (C . 120111 in prepara- 
tion) andlAkivam a et al. l (M. I201ll in preparation). In 
the ELAIS-S1 field we use the spectroscopic red shifts of 
X-ray sources presented bv lFeruglio et al.l (|2008ft . Over- 
all this provides high resolution spectroscopic redshifts 
for 270 of the X-ray sources within PRIMUS. Approxi- 
mately 50% of these sources are classified as BLAGNs. 
For the sources with galaxy classifications we determine 
the precision of the robust PRIMUS redshifts to be 
er 2 /(l + z) = 0.008. The outlier rate depends on mag- 
nitude: at i < 22, corresponding to 85% of our final 
sample (see Section @] below), the outlier rate (fraction 
with 6z/(l + z) > 0.1) is 4.8%; at fainter magnitudes the 
outlier rate may rise to ~ 30%, but is fairly poorly deter- 
mined due to the lack of sources reliable high-resolution 
spectroscopic follow-up to this depth. We note that the 
outlier rate for the X-ray detected galaxy population is 
slightly higher than determi ned for the over all PRIMUS 
galaxy sample (2.0%: see iCoil et all 12011ft but is low 
enough that we do not expect a significant impact on 
our results. 

2.2.1. XMM-LSS 

In t he XMM-LSS field we use the catalog of lUeda et all 
(2008) from the SXDS, based on seven deep XMM- 
Newton pointings (1.14 deg 2 ) that overlap our PRIMUS 
observations. We suppl ement this with t he point source 
catalog published in iPierre et al.l (|2007| ) , which utilizes 
the first 45 pointings (5.5 deg 2 ) of the XMM-LSS X-ray 
survey, although the overlap with PRIMUS is limited. 
Our combined catalog contains 4466 unique sources down 
to hard band (2-10 keV) fluxes of / 2 -iokcV ~2x 10~ 15 
erg s _1 cm -2 , 1073 of which fall within the area covered 
by PRIMUS masks, and 617 of which have secure opti- 
cal counterparts. We obtained robust redshifts and clas- 
sifications for 266 of these X-ray sources with PRIMUS 



over a total area of 1.78 deg 2 that is covered by both 
our PRIMUS masks and the X-ray data. We note that 
an additional ~ 1 deg 2 of our PRIMUS area is covered 
by X-ray data from the XMM-LSS survey that was not 
included in the IPierre et al.l ()2007ft catalog, and is thus 
not used here. 

2.2.2. COSMOS 

The entire 2 deg 2 COSMOS field has b een observed 
with XMM-Newton (jHasinger et al.l 12007ft to depths of 
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3 x 10 ergs cm , much deeper Chan 
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dra data reaching /2-iokcV ~ 8 x 10 
were also obtained for the central ~ 0.9 deg 2 ([E lvis ct al. 
2009ft . We merged t he XMM-Newton (jCappelluti et al.l 
2009) and Chandra (|Elvis et al.l [2009) point source cat- 
alogs, and removed duplicates by first matching bright 
sources (/0.5-iokoV > 10~ 14 ergs~ 1 cm -2 ) to correct 
for any overall offset between the two catalogs, and 
then searching for the counterparts of the XMM-Newton 
sources within the Chandra data using a 5" search radius. 
For matched sources we adopt the positions and fluxes 
from the Chandra catalog. We find that there are ~ 70 
XMM-Newton sources that fall within the approximate 
area observed with C handra but l ack C handra counter- 
parts, consistent with lElvis et al.l ()2009ft and attributed 
to areas or small gaps with low Chandra exposure, or 
some expected fraction of spurious sources. We retain 
all unmatched XMM-Newton sources in our final merged 
catalogs; spurious sources are very unlikely to have secure 
optical counterparts above our PRIMUS magnitude lim- 
its so will not significantly contaminate our sample. Our 
merged catalog contains 1232 unique hard X-ray sources 
over the 1.07 deg 2 with both PRIMUS and X-ray cover- 
age; 307 of these have robust optical counterparts, were 
targeted by PRIMUS, and have robust redshift measure- 
ments and classifications from the PRIMUS data. 

2.2.3. ELAIS-S1 



IPuccetti et al.1 (|2006ft presented a point source catalog 
from a mosaic of four partially overlapping XMM-Newton 
pointings in the ELAIS-S1 field that reach depths of 
/2-iokcV ~2x 10 -15 erg s _1 cm -2 . This provides X-ray 
data coverage for 0.52 deg 2 of our total 0.9 deg 2 PRIMUS 
area, although PRIMUS has little overlap with the deep- 
est portions of the ELAIS-S1 X-ray data. One hun- 
dred and twenty-nine hard X-ray sources fall within the 
PRIMUS area of which 105 have secure optical counter- 
parts. Robust redshifts and classifications for 59 sources 
are provided by PRIMUS. 

3. STELLAR MASS ESTIMATES 

We estimate stellar masses by fitting the ob- 
served optical/ncar-infrarcd spectral energy distributions 
(SEDs) of the galaxies in our sample. In the XMM- 
LSS field we construct observed SEDs based on the 
u*g'r'i' ' z' imaging from M egacam on the Cana da- France 
Hawaii Telescope (CFHT: iBoulade et"aTll2003ft obtained 
as part of the CFHTLS-Wide survejEH we use the pho- 
tometric catalo gs released by the C FHTLS Archive Re- 
search Survey (JErben et al.ll2009ft . In the COSMOS 



field we use BjVjg + r + i + z + imaging obtained with the 



1 http://www.nottingham.ac.uk/astronomy/UDS/data/data.html 
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TABLE 1 
Details and Source Numbers for our PRIMUS/X-ray fields 



Field Name 


Area (deg 2 ) 
(1) 


(2) 


X- 

(3) 


ray Source Numbers 
(4) (5) (6) 


(7) 


(8) 


Number of Galaxies 
(9) 


XMM-LSS 

COSMOS 

ELAIS-S1 


1.78 
1.07 
0.52 


1073 
1232 
129 


617 
914 
105 


410 
514 

82 


266 
307 
59 


116 

161 
25 


102 
111 
22 


83 
138 
21 


12523 
10114 
3079 


Total 


3.36 


2434 


1636 


1006 


632 


302 


268 


242 


25716 



Column 1: area with joint X-ray and PRIMUS coverage. 

Column 2: total number of hard-band detected X-ray sources in X-ray/PRIMUS area. 

Column 3: number with robust optical counterparts. 

Column 4: number of counterparts targeted by PRIMUS. 

Column 5: number with robust redshifts (Q > 3) from PRIMUS. 

Column 6: number with robust redshifts in the range 0.2 < z < 1.0 and X-ray luminosities 

in the range 42 < log Lx < 44. 

Column 7: number with robust redshifts, 0.2 < z < 1.0, 42 < logLx < 44 and galaxy 

classifications. 

Column 8: final number with robust redshifts, 0.2 < z < 1.0, 42 < logLx < 44, galaxy 

clas sifications and .M* above the redshift-dependent stellar mass limit determined in Section 

\n\ 

Column 9: total number of galaxies in sample with robust redshifts, 0.2 < z < 1.0 and M* 
above the redshift-dependent stellar mass limit. 



Suprime-Cam instrum ent on the Subaru 8.2m telescope 
(|Mivazaki et all [20021 ). u* and i* imaging from Mega- 
cam on the CFHT, and K s imaging from the Wide- 
field Infrared Camera at the CFHT (jMcCracken et alj 
120101 ): photometry in all these bands is obtained from 
the public multiwavelength photometric catalog released 
in 2009 ApriQ In the ELAIS-S1 field we use BVR 
imaging obtained with the Wide Field Imager at th e 2.2 
m La Silla ESO-MPI telescope (|Berta et all 2006) and 
Iz imaging based on observations with the Visible Multi 
Object Spectrogra ph camera at the Very Large Telescope 
(jBerta et al.ll2008f l. 

We perform the SED fitting using iSEDfit (see J. 
Moustakas et al., in preparation, for complete details). 
Briefly iSEDf it is a Bayesian SED-fitting code that com- 
pares the observed photometry against a large Monte 
Carlo suite of model SEDs spanning a wide range of age, 
metal licity, dust content, and star formation histor y (see, 
e.g., iKauffmann et al.l I2003b1 : iSalim et all 120071) . By 
marginalizing over all other parameters, we construct 
the posterior probability distribution of stellar masses, 
p(M*), for each galaxy encapsulating both the photo- 
metric uncertainties and the physical degeneracies among 
different models. We take the median of the p(A4*) dis- 
tribution as our best estimate of the stellar mass. 

We construct our models using the lBruzual fc CharlotJ 
(120031) stellar p opulation synthesis models, assuming the 
iChabrierl (|2003l ) initial mass function (IMF) from 0.1 to 
100 A^q. We assume a uniform prior on stellar metallic- 
ity in the range 0.004 < Z < 0.05, and smooth, exponen- 
tially declining star formation histories, ij){t) oc exp~ 7 *, 
with 7 drawn uniformly from the interval [0.01, 1] Gyr -1 . 
The time since the onset of star formation is given by t, 
which we draw from a uniform distribution ranging from 
0.1 to 13 Gyr. However, we apply an additional prior 
that disallows ages older than the age of the universe at 

1 http://irsa.ipac.caltech.edu/data/COSMOS 



the redshift of each galaxy. Final ly, we adopt the time- 
dependent dust attenuation law of lCharlot fc Fal] (2000) 
in which stars older than 10 Myr are attenuated by a fac- 
tor fi less than younger stars. We draw the U-band op- 
tical depths from an order two Gamma distribution that 
peaks around Ay « 1.2 mag, with a tail to Ay ss 6 mag, 
and n from an order four Gamma distribution that ranges 
fro m zero to unity centered on a typical va lue (/i) =0.3 
fsee lCharlot fc Falll[2000llWild et alJ[20Tll) . 

Adopting different prior distributions — for example, 
superposing stochastic bursts on otherwise smooth star 
formation histories, or adopting different metallicity 
or reddening priors — changes our stellar masses by < 
0.1 dex in the mean, and by typically less than a fac- 
tor of two for individual objects (see J. Moustakas et al., 
in preparation, for details) . Given that our sample spans 
more than two orders of magnitude in stellar mass, we 
do not expect these uncertainties to affect our conclu- 
sions. The largest potential systematic effect is due to 
our assumed IMF; however, that will simply shift all our 
stellar mass estimates by an equal amount. 
4. SAMPLE SELECTION 

In this paper, we wish to investigate the fraction of 
X-ray AGNs within the PRIMUS galaxy sample as a 
function of both the galaxy properties (A4*, color) and 
the X-ray luminosity. Our galaxy sample is restricted to 
(1) sources within the 3.37 deg 2 of PRIMUS with X-ray 
coverage; (2) ro bust redshift measurements (Q > 3, see 
iCoil et aLll2011l and R. J. Cool et al. in preparation for 
details); and (3) galaxy classifications - objects where 
the PRIMUS spectrum is best fit by a BLAGN template 
are excluded (see Section |4~31 below) . Our X-ray data are 
then used to identify AGNs within this galaxy sample. 

Our sample includes both primary PRIMUS tar- 
gets, with well defined targeting (sparse and density- 
dependent) weights, and a small number of secondary 
targets that are drawn from the same parent sample of 
sources but do not have a well-defined weight. We ex- 
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perimented with further restricting our sample to only 
primary targets, and applying corrections based on the 
targeting weights. This reduces our sample of X-ray de- 
tected AGNs but did not otherwise significantly alter the 
results presented in the rest of this paper. This indicates 
that the targeting incompleteness of the PRIMUS sample 
affects our galaxy sample and X-ray detected AGN sub- 
sample roughly equally; the small scale collisions that de- 
termine the targeting weights do not impact our derived 
fraction of AGNs. We therefore retain both primary and 
secondary targets to maximize our sample sizes. 

When a reliable high-resolution spectroscopic redshift 
is available from one of the catalogs described in Section 
12.21 above, we adopt this value rather than the PRIMUS 
redshift. In the COSMOS field we adopt reliable redshifts 
from the zCOSMOS-lOk bright sample (|Lilly et al.ll2009[ ) 
with a confidence flag of 3 or 4; we d o not use th e addi - 
tional redshifts of X-ray sources from lBrusa et al.l (|2010f ) 
in place of our PRIMUS redshifts due to the heterogene- 
ity of this catalog and the lack of the redshift confidence 
flags for individual sources. This provides updated red- 
shifts for 1643 sources in our three PRIMUS fields with 
X-ray coverage. 

TableQ]gives the numbers of galaxies and X-ray sources 
in our three fields that satisfy our various selection cuts. 

4.1. X-ray selection 

We restrict our AGN identifications to hard-band (2- 
10 keV) detections. Hard-band selection identifies a sam- 
ple that includes both X-ray-unabsorbed and moderately 
absorbed AGNs, although it will miss the most heav- 
ily absorbed, Compton-thick sources. Soft (0.5-2 keV) 
or full (0.5-10 keV) band selection, by comparison, will 
be biased toward the detection of unabsorbed sources 
(equivalent hydrogen column densities N-& ;$ 10 22 cm -2 ). 
Hard-band selection also ensures that we can accurately 
measure a rest-frame 2-10 keV X-ray luminosity based 
on the observed hard-band flux; we calculate luminosi- 
ties assuming a photon index r = 1.9. A disadvantage 
of hard band selection is the reduction in the size of our 
X-ray detected sample due to the lower sensitivity of X- 
ray observatories at higher energies. We further limit 
our AGN sample to Lx > 10 42 erg s _1 to ensure that 
our measured X-ray luminosity can be associated with 
an AGN (rather than star formation), and Lx < 10 44 
erg s _1 to exclude the brightest AGNs that may con- 
taminate the host galaxy light and to reduce the effects 
of our selection against unobscured AGNs (see Section 
14.31 below). Very few of the X-ray detected AGNs within 
our galaxy sample lie above Lx = 10 44 erg s _1 so this 
cut has little effect on our results. Hard band selection 
also ensures our AGN selection criteria are not indirectly 
affected by the host galaxy; conversely, if we choose to 
include soft-band-selected sources we could be biased by 
absorption effects within the host galaxy 

The X-ray selection of AGNs will be subject to signifi- 
cant incompleteness in the data that we must correct for. 
The depth of our X-ray data varies between our fields and 
within a single field due to the differing X-ray observa- 
tories, the varying vignetting and point spread functions 
within a pointing, and the different exposure times of the 
pointings. Our entire survey area is sensitive to bright 
fluxes, /2-iokcV ^2x 10~ 14 erg s _1 cm -2 , thus we will 
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Fig. 1. — X-ray area curves for each of our fields and total com- 
bined sample derived by comparing the number of sources in our 
compiled catalogs at a given hard (2-10 keV) X -ray flux with pre- 
dictions based on the log N — log S relations ofTGcorgakak is et all 
(2008). These allow us to correct for incompleteness due to the 
varying depths of our X-ray data, both between our fields and 
within an individual field. 

identify all AGNs with L x > 10 44 erg s" 1 out to z = 1 
(with the exception of Compton-thick sources) . However, 
this luminosity cut would severely reduce our sample size 
and we are mainly interested in studying the moderate 
luminosity (Lx ~ 10 42 ~ 44 erg s" 1 ) AGN pop ulation that 
domi nates the luminosity density at z < 1.0 ([Aird et al.l 
12010ft and where we are able to probe the properties of 
the host galaxy (as for the majority of this population 
the AGNs are moderately obscured). Instead of the strict 
flux cut, we apply completeness corrections to our de- 
rived AGN fractions based on the X-ray area (sensitivity) 
curves. We calculate the area curves by comparing the 
number of hard X-ray detected sources in bins of 0.1 dex 
from our parent point source catalogs (regardless of op- 
tical counterparts or PRIMUS targeting) that fall within 
our PRIMUS window with the predicted numbers as a 
fu nction of flux based on th e X-ray log N — log S relation 
of Georgakakis et al. l (|2008l ) . We fix the area curve at the 
total area coverage above the lowest flux where > 99% of 
the predicted sources are detected. This prevents small 
number statistics adding noise to our area curves at the 
highest fluxes. Area curves derived in this way are shown 
in Figure [T] 

Two different approaches to incorporate this complete- 
ness information are described in Sections [5] and [6] below. 
Figure [2] shows the distribution of X-ray luminosities for 
our X-ray detected sample, and the bins applied in Sec- 
tion[5j The effects of the X-ray incompleteness are clearly 
seen in the lack of lower Lx objects at high redshifts. 

4.2. Stellar mass limits 

The galaxy sample will be subject to a number of in- 
completenesses due to the targeting and redshift success 
rates that will depend on the optical apparent magni- 
tude and redshift, and thus indirectly depend on stellar 
mass. However, the AGNs will be subject to the same 
incompleteness effects (provided that the optical light is 
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Fig. 2. — X-ray luminosity (2-10 keV) vs. redshift for the X- 
ray-detected sub-sample of galaxies (crosses) and X-ray detected 
broad-line AGsN (red diamonds). Dashed lines indicate the bins 
used in Section [5] We restrict our sample to moderate-luminosity 
AGNs (42 < logLx < 44) below the break in the X-ray luminos- 
ity function. Our sample suffers from incompleteness due to the 
varying X-ray flux limits over our survey, which is particularly se- 
vere for logLx ;$ 42.5 at higher redshifts; this is fully described by 
our X-ray area curves (Figure [TJ and accounted for by our X-ray 
completeness corrections described in Sections \E\ and [6] We also 
exclude the broad-line AGNs, which are mainly found at higher 
X-ray luminosities (logLx 2^ 43.5). 

dominated by the host galaxy and not the AGN). Thus 
the fraction of AGNs within a sample of PRIMUS galax- 
ies for a limited stellar mass and redshift range should 
reflect the true fraction of all galaxies that host AGNs 
at that stellar mass and redshift. Nevertheless, we must 
set a lower limit on the stellar masses of galaxies in our 
sample to limit strong incompleteness effects that could 
depend on the host galaxy color and could bias our mea- 
surements of the AGN fraction. First, we classify the 
galaxies as blue or red if they lie above or below the 
following relation that defines the "green valley" , 



u-g = 0.671 - 0.031M, - 0.065z 



(1) 



where 



!l 



and M„ are the rest-frame colors 



and absolute magnitud es (derived using KCORRECT; 
iBlanton fc RoweisTl2007t ) and z is the redshift from 
PRIMUS. This cut is determined by identifying the peaks 
in u — g color that correspond to the red sequence in 
0.5 mag wide bins (in the range —21 < M g < —19) 
for galaxies with PRIMUS redshifts 0.4 < z < 0.7. A 
cut is made close to the minimum of the distributions 
and Gaussian fits are used to accurately determine the 
positions of these peaks. A linear fit is performed to de- 
termine the slope of the red sequence as a function of 
M g . We then take all galaxies in bins of Az = 0.1 over 
the range 0.4 < z < 0.7 and adjust this slope blueward 
until it passes through the minimum corresponding to 
the green valley. We fit this color offset as a function of 
redshift, which allows us to account for the overall color 
evolution of the galaxy population. Finally we identify 
the overall offset of the green valley for all galaxies at 
0.2 < z < 1.0 after accounting for the redshift evolution 
and magnitude dependence. 
In Figure [3] we show the distribution of stellar masses 
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Fig. 3. — Stellar mass vs. redshift for our PRIMUS galaxy 
sample, divided into red (top) and blue (bottom) galaxies accord- 
ing to Equation [T] Contour levels enclose 30%, 60% and 90% of 
the galaxy sample. X-ray-detected sources are shown by crosses. 
We define redshift-dependent stellar mass limits for the COSMOS 
(solid black curve), XMM-LSS (dashed black curve), and ELAIS- 
Sl (dotted black curve) fields based on the ma ss-to -light ratio of 
a maximally old stellar population (see Section I4.21 ) and the lim- 
iting targeting magnitude in each field. We cut both the red and 
blue galaxy samples using these mass limits; for blue galaxies we 
are able to probe to lower masses but we wish to limit any color- 
dependent incompleteness effects and thus apply this stricter cut. 
The dashed gray lines indicate the bins used in Section [5] 

as a function of redshift for our red and blue galaxy 
samples and the X-ray detected sub-samples. We de- 
fine a redshift-dependent stellar mass limit for each field 
above which strong, color-dependent incompleteness ef- 
fects should be minimized. To define this limit we take 
a maximally old (i.e., r = 0) simple stellar population 
with a formation redshift of z = 5 and determine mass- 
to-light ratios (for the targeting optical bandpass) as a 
function of redshift. We use these mass-to-light ratios to 
convert our targeting magnitude limits in each field to 
an effective stellar mass limit as a function of redshift. 
These stellar mass limits are shown by the black curves 
in Figure G3 Above this mass limit our sample should 
be representative of the entire galaxy population as no 
stellar population should be redder than our maximally 
old template. We apply this cut to both the red and blue 
galaxy samples. We are able to detect blue galaxies and 
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indeed some red galaxies down to much lower masses, 
but above this mass cut we will detect both red and blue 
galaxies. This ensures that we can accurately measure 
the probability of a galaxy of a given stellar mass hosting 
an AGN, and will not be biased if AGNs preferentially 
reside in red or blue host galaxies. In section 18.11 we 
will investigate the positions of AGN host galaxies in the 
color-magnitude diagram in greater detail. 

The dashed grey lines in Figure [3] show the redshift 
and mass bins used to derive the AGN fraction in sec- 
tion [5] The mass cut passes through our lowest stellar 
mass bins, reducing the number of galaxies in that bin. 
The number of detected AGNs should also be reduced by 
the same amount, and thus we will accurately measure 
the AGN fraction, assuming the AGN fraction does not 
strongly depend on mass or redshift within that bin. In 
Section |6] we use a maximum likelihood approach to fit 
functional forms to the AGN fraction. The same mass 
cut is applied, but no binning is required and thus any 
biases associated with the assumption of a constant frac- 
tion within a bin are avoided. 

4.3. Broad-line AGNs 

Our final selection criterion requires a galaxy classi- 
fication based on the PRIMUS template fit. This will 
exclude unobscured AGNs that exhibit broad emission 
lines in the optical spectrum. We classify objects based 
on the relative \ 2 values for the best-fit BLAGN and 
galaxy templates (see R. J. Cool et al., in preparation, 
for details). This ensures we only exclude objects that 
are significantly better fit by a BLAGN template. 

The BLAGN cut is essential as we can only accurately 
measure stellar masses for AGNs where the host galaxy 
dominates the optical light. This does introduce an addi- 
tional incompleteness that we do not attempt to correct, 
and thus our study of the AGN fraction is only applica- 
ble to the "obscured" population. Our definition of "ob- 
scured" means simply lacking strong broad lines in the 
optical spectrum; this could be due to intrinsic obscura- 
tion close to the SMBH (the torus), obscuration within 
the host galaxy or dilution of the light from a weak 
AGN within a substantially brighter host galaxy. Fig- 
ure [5] shows that the BLAGNs are mainly found at high 
X-ray luminosities, Lx > 10 435 erg s _1 , consistent with 
the increase in the fraction of uno bscured AGN s with 
increasing X-ray luminosity (e.g., iTreister et al.l 120091 : 
iTrump et al.ll2009h . We restrict our sample to Lx < 10 44 
erg s _1 to limit the effects of this incompleteness, but we 
do not reduce the limit further as this would decrease 
our sample size and dynamic range. 

5. THE PREVALENCE AND DISTRIBUTION OF AGN 
ACCRETION ACTIVITY AS A FUNCTION OF STELLAR 

MASS 

In this section we investigate how the prevalence and 
distribution of AGN accretion activity, traced by the 
number of X-ray-detected AGNs and their X-ray lu- 
minosities, depends on the stellar masses of potential 
host galaxies. We consider two broad redshift bins that 
evenly divide the redshift range probed by PRIMUS, 
0.2 < z < 0.6 and 0.6 < z < 1.0 as shown in Figures 
[5] and |3J Our choice of two redshift bins is a compro- 
mise between including enough objects to reveal trends 
with luminosity or mass, yet avoid serious biases due to 



possible redshift-dependent completeness or evolutionary 
effects (evolution of the fraction with redshift is investi- 
gated in Section [6]). 

Within each redshift bin we subdivide our galaxy sam- 
ple into bins of different A4*; we then identify X-ray 
AGNs within each of these stellar mass bins. The dis- 
tribution of X-ray luminosities, Lx, for objects within 
this A4* bin is given by p(Lx \ A4*,z), the conditional 
probability density function describing the probability of 
a galaxy of a given .M* and z hosting an AGN with 
luminosity Lx- We define p(Lx \ A4*,z) as the probabil- 
ity density per logarithmic luminosity interval, therefore 
with units of dex" 1 . We use the term "the prevalence 
of AGN activity" to refer to the probability of a galaxy 
hosting an AGN defined by this function, whereas "the 
distribution of AGN activity" refers to the distribution of 
black hole growth traced by the luminosity dependence 
of this function. The goal of this section is to measure 
p(Lx I M*, z) in our two redshift bins. 

We can relate p(Lx \ M.*,z) to the AGN fraction, 
/agn, which we define as the fraction of galaxies of a 
given stellar mass and at a given redshift that host an 
X-ray (non-broad- line) AGN with Lx > 10 42 ergs -1 , 
thus 



/agn(A4*,z) 



p(Lx\M, 1 z)dlogLx. (2) 



;L X =42 



We note that a basic "AGN fraction", the fraction of 
galaxies with an AGN, is a poorly defined and somewhat 
ad hoc quantity that can depend on the depths of the 
data and the method used to identify AGNs, in addi- 
tion to the selection of the parent galaxy sample. This is 
certainly true when using X-ray data to identify AGNs; 
the fraction of galaxies that are X-ray sources is highly 
dependent on the depths of the X-ray data and the red- 
shifts and magnitude limits of the parent galaxy sample. 
Our quantity of interest, p(Lx \ M*, z), is clearly defined 
and contains additional information on the distribution 
as a function of X-ray luminosity. 

5.1. Completeness weighting 

For our initial investigations we bin our sample by red- 
shift, stellar mass, and X-ray luminosity, and directly 
measure p{Lx \ M.*,z) in these bins. The observed dis- 
tribution of X-ray luminosities for a given M.* and z 
bin will not directly correspond to p{Lx \ M*,z) as the 
varying sensitivity of our X-ray data prevents us from 
detecting low-luminosity or high-redshift X-ray AGNs 
in all galaxies of a given A4* within our sample. It is 
thus vital to correct for these effects to accurately mea- 
sure p(Lx I M*,z) and reveal the relationship between 
AGN activity and the stellar masses of their host galax- 
ies. To perform X-ray completeness corrections within 
each bin we first determine the total number of galax- 
ies, iVg a i, within a mass and redshift bin shown in Figure 
[3) We then identify X-ray-detected AGNs within this 
same mass and redshift bin. We take the X-ray luminos- 
ity of each source, L^, and calculate Pdct{Lk , Zj) — the 
probability of being able to detect an AGN with this 
X-ray luminosity Lk at the redshift Zj of each galaxy if 
the galaxy hosted an AGN of this luminosity. This is 
determined by the sensitivity of our X-ray data. The 
area of our survey sensitive to a given flux, A(fx), is de- 
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scribed by the X-ray area curves determined in Section 
Ol Thus, 



Pdet(L k , Zj) 



A{f x (L k , Zj )) 



itotal 



(3) 



where -A to tai = 3.36 deg 2 , the total area covered by both 
PRIMUS and our X-ray observations, and fx{L k ,Zj) is 
the hard band X-ray flux that would be observed from 
a source of luminosity L k at redshift Zj . We sum these 
probabilities over all galaxies in our mass and redshift 
bin to determine a weight for each X-ray source, 



N, 



w k 



al 



•A', 



J2j=f Pdet(L k , Zj) 



(4) 



where w k is the weight for the fcth X-ray-detected AGN. 
Finally, we obtain completeness-corrected estimates of 
p(Lx I M.*,z) in a range of X-ray luminosity bins by 
summing these weights within each luminosity bin and 
dividing by the total number of galaxies (in our given 
stellar mass and redshift bin). Thus, 



A. 



p{L m I M*,z) 



J2 Wk 

fc=i 



A ga i A log L r , 



(5) 



where we sum over all A,^ GN AGNs in the mth X-ray lu- 
minosity bin of width A log L m . We calculate an error on 
p{L m J M r , z) based on the Poisson error in the number 
of X-ray detected sources in a bin (applying the formulae 
of iGeh rcls 1986J); we do not consider the contribution of 
the statistical error in the total number of galaxies in 
a bin which is minimal due to the very large numbers 
of sources compared to the X-ray-detected population. 
Our completeness corrections make the reasonable as- 
sumption that the AGNs (and their host galaxies) that 
are not X-ray detected due to incompleteness have the 
same properties as the sources that are detected at the 
same Lx, M.*, and z. 

5.2. Results 

Our estimates of p(Lx \ M*,z) for our two redshift 
bins are shown in Figure |4] where different colors indicate 
different stellar mass bins and different symbols indicate 
different X-ray luminosity bins. Two general trends im- 
mediately apparent in both our redshift bins are as fol- 
lows: 

1. The probability of a galaxy hosting an AGN of 
a given X-ray luminosity increases significantly as 
stellar mass increases. 

2. The probability of a galaxy of a given stellar mass 
hosting an AGN decreases as X-ray luminosity in- 
creases. 

A number of studies have shown that AGNs are found 
predominately in massive galaxies (e .g. JBest et al.ll2005t 
IBrusaet alj|2009al IXue et alj|2010h . We also demon- 
strate that the probability of hosting an AGN de- 
creases with increasing X-ray luminosity at all stel- 
lar masses. Luminous AGNs are much rarer than 
their low-to-moderate-luminosity equivalents, as shown 



by t he AGN luminosity function for th e overall popula- 
tion (jUeda et alJl2003tlAird et al.ll2010j) . thus the overall 
trend of fewer luminous AGNs is expected. However, 
more luminous AGNs may be expected to be associated 
with the most massive galaxies, while low-mass galaxies 
might be expected to host the low-luminosity AGN pop- 
ulation. Our results show that this is not the case. In- 
deed, for a sample of galaxies of any given mass, a higher 
fraction will be found to host low-luminosity AGNs. The 
X-ray luminosity distribution of AGNs does not appear 
to depend strongly on the host galaxy mass; equivalently 
the probability of hosting an AGN as a function of host 
galaxy mass appears to be independent of the X-ray lu- 
minosity. 

To further investigate these trends we perform \ 2 fits 
to our AGN fractions using our binned estimates and 
Poisson errors. Motivated by the distributions seen 
in Figure HI which show roughly linear trends in log- 
log space, we assume simple power-law relations for 
p(Lx I M.*,z) as a function of stellar mass at a fixed 
X-ray luminosity, 



log [p{L x I M*,z)] = a + 6 log 



M> 



10 n Af m 



(6) 



where a and b are the intercept and power-law slope of 
the relation. The best-fit values for each Lx bin are 
shown in the small insets on the left side of Figure |H 
The fits are consistent with a slope that does not change 
with X-ray luminosity, in agreement with our basic obser- 
vation that the variation of p(Lx \ M.*, z) as a function 
of M.* is independent of Lx- We also fit power-law rela- 
tions for p{Lx I M.*,z) as a function of X-ray luminosity 
for a fixed stellar mass, 



\og\p{Lx\M*,z)] 



dlog 



Lx 



lO^ergs- 1 



(7) 



with intercept c and slope d. The best fit parameters for 
each stellar mass bin are shown in the inset panels on the 
right side of Figure [4] Again, the fits are consistent with 
the distribution of AGN luminosities being independent 
of stellar mass. 

We note that the same trends are seen in both redshift 
bins. Given large errors, the slopes of our fitted power 
laws are consistent between our redshift bins. There is 
also weak evidence that the intercepts, shown in the inset 
panels, are generally larger at higher redshift, indicating 
evolution with redshift in the overall fraction. This is 
investigated in more detail in the next section. 

6. REDSHIFT EVOLUTION 

To investigate any evolution oip(Lx \ M.*, z) with red- 
shift we must divide our sample into finer redshift bins. 
This reduces the number of sources in each of our stel- 
lar mass and luminosity bins. The \ 2 fitting approach 
requires fairly heavy binning and will be sensitive to the 
chosen bin sizes. If the bins are too large then the actual 
trends in the data may be hidden, and the individual 
binned measurements may be systematically biased by 
large sensitivity variations within a bin; if the bins are 
too small then they will contain too few objects and x 2 
statistics cannot be used. Binning also means that we do 
not utilize all the available information: the actual mea- 
surements of A4* and redshift for each individual galaxy 
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Fig. 4. — Estimates of p(Lx I M*,z), the probability density for a galaxy of given stellar mass, Al*, and redshift, z, to host an AGN of 
X-ray luminosity, Lx- Integrating this function over Lx gives the AGN fraction at a given stellar mass and redshift (see Section [5] Equation 
[2t . Results are shown as a function of Al* (left panels) and Lx (right panels) for two redshift ranges: 0.2 < z < 0.6 (top) and 0.6 < z < 1.0 
(bottom). Colors/symbol sizes indicate the Al* bin and symbol types indicate the Lx bin in both the left and right panels. Error bars are 
based on the Poisson error in the number of X-ray detected AGN within the bin. Lines/small symbols show best fit power-law relations for 
p(L^ | Al,, z)as a function of Al* (left) or Lx(right) for each Lx bin (left, different symbols) or A4» bin (right, different colors) as defined 
in Equations ]E\ and [7] respectively. The inset panels show the best-fit slopes and intercepts of these power-law relations. The horizontal 
black dashed lines show the weighted means of the slopes in the different Al* and Lx bins. We find that p(L^ \ Al*, z) increases at higher 
Al*, but the slope of the relation is consistent in each Lx bin (left panels); p(Lx | M*, z) also falls with increasing Lx, but the slope does 
not depend on the Al* bin (right panels). 



and the Lx measurements for each X-ray detected AGN. 
Wc therefore develop a maximum likelihood fitting ap- 
proach to determine the dependence of p(L x \ M.*-, z) on 
A/f* and Lx in a number of smaller redshift bins. In Sec- 
tion 16.31 we directly incorporate redshift evolution into 
this scheme, removing the need to bin by redshift. 

6.1. Maximum likelihood fitting 

The results of the previous section imply that 
p(Lx | M*,z) at a fixed redshift may be written as a 
separable function of stellar mass and X-ray luminosity 
of the form 

p(L x \M*,z)d\ogL x = K(j^) f-^J dlogLx, 

where A4q and Lq are arbitrary scaling factors, K is 
an overall normalization and jm and jl are power-law 



slopes. We adopt this functional form to perform an un- 
binncd extended maximum-likelihood fit to our data, and 
measure values of our three free parameters (K, jm and 
7l) in each of our redshift bins. The extended maximum- 
likelihood fit allows us to directly constrain the normal- 
ization, K, along with the other model parameters (jl, 
^m) and estimate robust confidence intervals on all three 
parameters. To find the best values we maximize the log- 
likelihood function, 



AGN 



ln£ 



-AT- 



N. 

E 



lnpfe 



(9) 



where the summation is performed over all JV^ X- 
ray-detected AGNs in the zth redshift bin, and pk is the 
probability that a galaxy of stellar mass Aik will host an 
AGN of X-ray luminosity L^ . This is directly related to 



PRIMUS: The dependence of AGN accretion on host stellar mass and color 



11 



M k V M (L k ^ L 



p(Lx | M*, z), thus 



Pk =p(L k \M k ,z k ) = K^ J fJ \^j , (10) 

where we assume no dependence on z within our fine 
rcdshift bin. The expected number of X-ray AGNs, N , 
for a given set of model parameters is found by sum- 
ming the probability of observing an AGN of any lu- 
minosity logLx > 42 over all Nf a galaxies in our red- 
shift bin. The probability of observing an AGN is the 
combination of the intrinsic probability density function, 
p(Lx | M*,z), which gives the probability that a galaxy 
hosts an AGN of luminosity Lx, and the probability of 
actually being able to detect such a source, Pdet(Lx,z), 
defined in Equation [3l Thus, 



N sal 
Af=V / p(L X | Xj,«j)Pdet(ix,2j)dl0g£ X (11) 



i=i 

jy-gal 



ma 1m r (l-^ il 



i_ 

Mo 



L 



Pdet(L x ,Zj) dlogix- 



The X-ray completeness is accounted for by the 
Pdet(Lx, Zj) term; this reduces the expected number of 
AGN at low luminosities and high redshifts due to the 
X-ray sensitivity limits of our survey. We set log L mm — 
42.0 as the cutoff luminosity above which an X-ray source 
is associated with AGN activity. We estimate equivalent 
la errors on the model parameters from the maximum 
projection of the AS = 1.0 likelihood surface for each 
parameter, where S = — 21n£. 



6.2. Results in finer redshift bins 

We have used the maximum likelihood fitting ap- 
proach described above to constrain the functional form 
of p{Lx | M*,z) in four redshift bins, evenly spaced 
across the 0.2 < z < 1.0 range probed by our sample. We 
fix the scaling factors at Mo = 1O 11 M and L a = 10 43 
erg s _1 . We note that our current scheme makes the sim- 
plifying assumption that p(Lx \ M*, z) does not strongly 
depend on redshift within the range of the redshift bin. 
The best fit parameters are shown by the points and er- 
ror bars in Figure [5] We find that the power-law slopes, 
jM and jl, do not change significantly with redshift. We 
do, however, find that the normalization increases signif- 
icantly with increasing redshift. Thus, the fraction of 
galaxies of a given stellar mass hosting AGN of a given 
X-ray luminosity increases substantially with redshift, 
indicating a higher prevalence of AGN accretion activity 
at earlier times. 



6.3. Incorporating redshift evolution 

We modify our maximum-likelihood estimator to in- 
clude redshift evolution, allowing us to directly constrain 
the form of the observed redshift evolution and remove 
redshift binning from our analysis. We assume that the 
redshift evolution is independent of X-ray luminosity or 
stellar mass and has a simple power-law form in (1 + z), 
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Fig. 5. — Best-fit parameters from the maximum-likelihood fit- 
ting for p(Lx I M.*i z). Squares show the best fit parameters (and 
la equivalent errors) fitting in four redshift bins, assuming that 
p(Lx | M*, z) is constant over the redshift bin. The dashed lines 
are the results fitting over the entire redshi ft r ange, allowing for 
redshift evolution as described by Equation 1121 the dashed lines 
in the top two panels show the best estimates of 7m and -y^ for 
the entire redshift range (gray regions indicate ltr uncertainty); in 
the bottom panel the dashed line indicates the e qui valent of K for 
the model with redshift evolution (see Equation ll3ll , We find that 
p(Lx | M*,z) is described by independent power-law functions of 
Mt and Lx that do not change with redshift; there is however 
significant evolution in the normalization that follows a power-law 
dependence on (1 + z). 



and thus we modify Equation [8] to 



p(L x | M*,z) dlogix = 
Mo) Uo 



(12) 



.4 



1 + 2 



1 



Z{) 



d log Lx 



where we set the scaling factor zq = 0.6, and A is 
an overall normalization. We adapt Equations [10] and 
1121 to incorporate this functional form and perform our 
maximum- likelihood fitting over our entire 0.2 < z < 1.0 
redshift range. Our best-fit parameters and their uncer- 
tainties are given in Table [21 We also show the best-fit 
values and uncertainties of the 7^1 and 7^ slopes in the 
top two panels of Figure [5] by the dashed lines and gray 
shaded regions respectively. In the bottom panel of Fig- 
ure [5] we show the equivalent normalization, 



K = A[{l + z)/{l + zo)V% 



(13) 
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TABLE 2 

Best-fit Parameters from Maximum-likelihood Fitting 

with Redshift Evolution 



Parameter 



Value 



logic A 

1M 



-1.75 
0.64 

-0.83 
3.67 



+0.03 

0.03 
+0.06 

0.05 
+0.06 

0.07 
+0.60 

0.61 



based on our best fit evolutionary model. We find evi- 
dence for very strong evolution in the AGN fraction with 
redshift of the form p(Lx | M*, z) oc (1 + z) 4 . 

To aid visualization of our data and comparisons 
with our models w e implement the N b s /N m< ±i method 
(|Mivaii et alj 12001). This method accounts for signifi- 
cant variations of both our model function and the X- 
ray sensitivity over a bin. This is especially useful given 
our relatively small sample of X-ray-detected AGNs, the 
large variation of the X-ray sensitivity as a function of 
Lx and redshift, and the strong redshift dependence we 
have found. This allows us to accurately combine data 
over large bins and identify true deviations of the data 
from the model. 

The Nobs/N m< n method compares the observed num- 
ber of sources in a bin, Nobs, with the predicted number 
based on a model fit. The binned estimate is then given 
by the model evaluated at the center of the bin, scaled 
by the ratio of N b s /N mi i\. For our purposes, JV m< u is 
the predicted number of galaxies within a given redshift 
and mass bin that we expect to detect as an X-ray AGN 
over a given range of luminosities, for a particular set of 
model parameters. Thus, 
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(14) 

where the summation is over all Nf a galaxies in the ith 
redshift and stellar mass bin, and £i£j and L 1 ^ are the 
upper and lower limits of the mth X-ray luminosity bin. 
In Figure [S] we show p(Lx | M.*,z) as a function of 
stellar mass and luminosity over the full redshift range 
(0.2 < z < 1.0) with Aobs/A^mdi binned estimates. The 
model lines are evaluated at z — 0.6, but our binned 
estimates use data over the full redshift range and are 
corrected for the strong redshift evolution we find in the 
AGN fraction. The binned estimates show that our ob- 
served data are well described by our model. 

7. EDDINGTON RATIO DISTRIBUTION 

In this section we investigate how the fraction of galax- 
ies hosting an obscured X-ray AGN depends on the Ed- 
dington ratio: the ratio of the bolometric luminosity 
to the Eddington limit, LEdd, where radiation pressure 
balances the infall of material under gravity (assuming 
spherical symmetry and hydrostatic equilibrium). The 



Eddington ratio tracks the rate at which matter is ac- 
creted by the SMBH and thus relates to the mode of 
accretion and the processes that fuel the AGN (e.g., 
IKauffmann fc Heckmanl [20091: iTrump et all l20llT) . The 
Eddington ratio could be the principal factor that deter- 
mines the prevalence of AGN activity in potential host 
galaxies. It is thus vital to accurately measure the dis- 
tribution of Eddington ratios, account for possible selec- 
tion effects that will bias against sources accreting at the 
lowest Eddington ratios, and uncover the relationship be- 
tween AGN activity and the properties and evolution of 
their host galaxies. 

7.1. Indications of an Eddington ratio dependence 

The results of Section [6] indicate that the prevalence 
and distribution of AGN activity is related to the Ed- 
dington ratio. We find that p(Lx | M*,z) increases 
with Al» at a fixed X-ray luminosity. However, the 
mass of the central SMBH that powers an AGN is cor- 
related with the mass of a galaxy's central stellar bulge, 
as shown by the Mhh — a relation (|Ferrarese fc Merritt] 
2000; Gcbhardt et al. 2000) as well as more direct probes 
(e.g., Hari ng fc Rixll2004[ ). A higher stellar mass galaxy 
is therefore expected to host a higher mass SMBH, 
and thus the same X-ray luminosity corresponds to a 
much lower Eddington ratio. Therefore, the rise of 
p(Lx | M*,z) with M* may simply reflect the higher 
luminosities of low Eddington ratio sources in galaxies 
with higher stellar masses. We also find roughly similar 
slopes for the dependence of p(Lx | M*,z) on M* and 
Lx such that 7l ~ —"/m to within ~ 2er. This indicates 
that p{Lx | M*, z) may be dependent on Eddington ra- 
tio, AEdd, and redshift only, assuming a single bolometric 
correction such that Lx°c Lboi and a direct proportion- 
ality between the stellar mass of a galaxy and its central 
SMBH, M*<x Mhh- Thus we can re-write Equation IT21 
as 

1 + zV* 



p(A E dd \M*,z) dlogA E dd = A Ag? ld ( j—— ) dlogA E dd, 



where je ~ 1l ~ -Jm ~ -0.7. 



(15) 



7.2. Bolometric corrections 



The assumption of a single bolo metric correcti o n may 
be a n over-simplifi c ation. Both IMarconi et al.1 (|2004T ) 
and [Hopkins e t al.1 (|200 7b) have shown that the 2— 
10 keV X-ray bolometric correction, fc2-iokcV, is lu- 
minosity dependent with an approximately power-law 
form, based on composite SEDs and the observed dis- 
tribution of chqx, the X-ray to optical spectral index 
(jZamorani et alll981f) . A number of recent studies, how- 
ever, have used observed multiband SEDs of unobscured 
AGNs to accurately measure bolometric luminosities and 
propose that fe-iokcV m &y be more closely correlated 
wit h the Eddington ratio, albeit with significant scat - 
ter (|Vasudevan L Fabian|[2007l l2009t Eusso et alJl20Toft . 
These studies focused on relatively luminous, unobscured 
AGNs, that are mostly accreting at high Eddington ra- 
tios (AEdd ^ 0.1), and the correlation docs not appear to 
hold at lower AEdd- We t herefore adopt t he lumi nosity- 
dependent /c 2 -iokcV from Hop kins et al.l (|2007bD . Our 
bolometric corrections vary from ~ 15 — 50, spannin 
the range of observed values seen in lVasudevan &; Fabia 
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Fig. 6. — Estimates of p(Lx I M*,z) as a function of A4* (left) and Lx (right) based on our maximum-likelihood fitting results. We 
combine data over our entire redshift range (0.2 < z < 1.0). The sm all g ray symbols and lines (left) and colored lines (right) show our best 
fit model from the unbinned maximum-likelihood fitting (Equation 1121 Table [2]l evaluated at the center of our redshift range (z = 0.6). 
Binned data points are estimated using the iVobs / ^mdl method (Section I6,3| l. which compares the observed number of AGNs with that 
predicted within an A4* and Lx bin using our best-fit model, given our observed parent galaxy sample. This allows us to correct for strong 
redshift, .M* and Lx dependent effects on the binned estimates, and directly compare deviations of the data from the model. Colors/symbol 
sizes and symbol types correspond to the M t and Lx bins, respectively (as in Figure [4}. Error bars are based on the Poisson error in the 
number of X-ray detected AGNs within the bin. 

pendence only are given in Table |31 In Figure [71 we 
show p(AEdd | M.*,z) as a function of AEdd for galaxies 
in different stellar mass bins; we combine our entire red- 
shift range and use our No^s/N^i approach described in 
Section HOI to account for the strong redshift evolution. 
Our results show that the distribution of Eddington ra- 
tios does not depend on the stellar mass of the galaxies; 
all of the points are consistent with the same global re- 
lation. For a sample of galaxies of any given mass we 
expect to find the same distribution of Eddington ratios. 
In the final panel of Figure UJ we combine data over our 
full range of stellar masses and redshifts, showing the ex- 
cellent agreement with our best-fit relation for the AEdd 
dependence. For comparison, we also show the observed 
distribution of Eddington ratios for our hard X-ray AGN 
sample without accounting for the X-ray sensitivity ef- 
fects or redshift evolution of p(AEdd | M*,z) (gray his- 
togram). This has a roughly lognormal shape that peaks 
at Asdd ~ 0.01 rather than the power-law shape that we 
recover. In Figure [8] we show p(AEdd | M.*,z) for our 
full stellar mass range but in different redshift bins; we 
find a significant increase in the fraction of galaxies that 
host AGNs with redshift, although the distribution of 
Eddington ratios follows the same power-law form. 

Our results indicate that AGNs are not predominantly 
in massive galaxies when considered in terms of Edding- 
ton ratio. Above a given X-ray luminosity, a higher frac- 
tion of AGNs will be found in more massive galaxies, but 
this is a selection effect as low Eddington ratio AGNs in 
high stellar mass galaxies will have higher X-ray lumi- 
nosities. We find that the Eddington ratio distribution 
is a universal function with a power-law shape of slope 
-0.65. 



(2007) for AEdd ^ 0.1, We do not expect uncertainty in 
the bolometric correction to significantly affect our re- 
sults given the wide range of Eddington ratios that we 
probe (~four orders of magnitude). 

We replace Lx with the bolometric luminosity, Lboi, in 
Equation 1121 and then repeat our maximum-likelihood 
fitting for the entire galaxy sample. Using this bolomet- 
ric correction we find a slightly flatter slope for the de- 
pendence on (bolometric) luminosity, jl — —0.70 ±0.05, 
which is in good agreement with the inverse of the mass- 
dependent slope, 7x = 0.64 ± 0.05, consistent with the 
prevalence of AGN activity being primarily determined 
by the Eddington ratio. 

7.3. Maximum-likelihood fitting of a universal 
Eddington ratio distribution 

We repeat our maximum-likelihood fitting for our en- 
tire galaxy sample using the form given in Equation [15l 
directly fitting a single power- law slope, 7^, for the Ed- 
dington ratio. We assume a constant scaling between 
the black hole m ass and the host galax y stellar mass, 
M hh w 0.002A4* ([Marconi fc Huntll2003l where we have 
assumed M* ~ .Mbuigc)- We use this relation to calcu- 
late Eddington ratios, 



AEdd — 



L\y 



ol 



^Edd 



(16) 



where the Eddington limit, L^dd, in units of erg s 1 
given by 



L R dd = 1.3xlO d 



M 



bh 



1.3 X 10 



:ix 



M Q 
0.002A/f* 



M 







(17) 



The results of our maximum-likelihood fitting for our 
entire sample with an Eddington ratio and redshift de- 



7.4. Specific accretion rates 

The Eddington ratio distributions presented in this 
section assume that the stellar mass of a galaxy is directly 
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or the redshift evolution of p(Agdd I M*,z). The observed distribution has a roughly lognormal shape with a peak at Asdd ~ 0.01 rather 
than the true power-law shape that we recover. 



proportional to the mass of the central SMBH. However, 
the A4bh — o~ relation implies a correlation between the 
mass of the stellar bulge and the SMBH. The fraction 
of the total stellar mass contained in the bulge will de- 
pend on the galaxy morphology, which in turn is related 
to both the color and s tellar mass of the host galaxy 
(see lBlanton fe Moi istakas 2009), introducing significant 
scatter to our A/fbh estimates. It is also unclear if the 
■Mbh — cr relation evolves with redshi ft or exists at all out- 
side the local universe (e.g., Shields et al. 2003; Merlonil 
2001 iCrotonl [20061: IPeng et alJl2006al IWoo et all 120061: 
Salviander et al.l 120071: IW00 et all 120081: iMerloni et all 
20101 : IShen fc Kellvll201fft . In fact. iJahnke et all (I2009f ) 
propose that while the .A/fbh — -Mbuige relationship may 
evolve with redshift, there is little evolution in the cor- 
relation between .A/fbh and the total stellar mass, Ai*. 
Nevertheless, uncertainties in our estimates of .A/fbh could 
affect our measured Eddington ratios and may bias our 
measurement of je or indeed our conclusion that the 
prevalence of AGNs is independent of the stellar mass 
of the host galaxy and primarily determined by the Ed- 
dington ratio distribution. It is perhaps surprising given 
these large uncertainties and systematic biases in our es- 
timates of Eddington ratios that we do recover a univer- 
sal Eddington ratio distribution with a constant slope 
that accurately describes our data throughout our red- 
shift range and over several orders of magnitude in .A/f* 



and A E dd- 

However, we can disregard any uncertainties in our 
inferred black hole masses if we consider our AEdd val- 
ues as a tracer of the specific accretion rate — the rate 
of black hole growth relative to the stellar mass of the 
host galaxy — rather than the true Eddington ratio. The 
results shown in Figure [7] show that the probability of 
finding an AGN with a specific accretion rate is described 
by a power-law distribution that is independent of stel- 
lar mass. Thus AGNs are not predominantly found in 
massive galaxies when we consider AGNs growing above 
a specific accretion rate. The conclusions laid out above 
and in Section [9] below are therefore robust despite the 
significant uncertainties in our estimates of A^bh- 

We note that evolution of the A4bh — o relation could 
alter the slope of the redshift e volution, •%- , that w e find 
for p(A Ed d | M^z). Indeed, IPeng et all (|2006bf ) find 
that black hole masses at z ~ 1 are a factor ~ 2 larger 
than implied by the local A^bh — o~ relation. However, 
significantly stronger evolution of the A/fbh — cr relation 
is required to account for the entire factor > 10 evolution 
we find in p(AEdd | A4*, z) between z ~ 0.2 and z ~ 1. 

In the remainder of this paper we discuss results in 
terms of Eddington ratio, but note that we actually mea- 
sure the specific accretion rate (with a particular scaling 
factor) and thus our results could instead be interpreted 
in terms of this quantity, providing more freedom in the 
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Fig. 8. — Best-fit model of p(Agdd I M*,z) evaluated at the 
centers of the four indicated redshift bins (lines) and compared to 
the data over the given redshift bin and our full range of stellar 
masses, 9.5 < log .M*/.M0 < 12.0 using the Afobs/^mdl method 
(points and error bars, slightly offset from the center of the A^dd 
bin to aid comparison). We find significant evidence for a strong 
evolution with redshift of the form p(Agdd I M*,z) oc (1 + z) 3,8 
(given our assumptions for the bolometric luminosity corrections 
and scaling of M t to A^bh)- However, with our data we cannot 
distinguish between an overall increase in the density of black holes 
that are actively accreting at a given A^dd with increasing redshift, 
or a shift of the e ntire population to higher Agdd ( see further dis- 
cussion in Section 19.31 1. 



TABLE 3 

Best-fit Parameters from Maximum-likelihood Fitting 

with Redshift Evolution and Eddington Ratio Dependence 



Parameter 



Value 



li ik 



IE 



n icr + 0.08 
■ J - lo -0.08 

-+0.04 
-0.04 



-0.65 1 



3.47 



+ 0.49 
-0.48 



underlying distribution of black hole masses. 

8. THE DEPENDENCE OF AGN ACTIVITY ON HOST 
GALAXY COLOR 

As discussed in the introduction, AGN accretion activ- 
ity and star formation processes appear to be connected. 
Both the overall star formation rate density and the AGN 
luminosity density have declined rapidly since z ~ 1, 
indicating that these processes may be triggered by a 
common mechanism. The strong correlation between 
black hole mass and the masses of central stellar bulges 
(trac e d by the Mhh — g r e lation: iFerrarese fc Merrittl 
120001 : iGebhardt etatl 120001: iTremaine et all I2002T ) also 
implies that these processes are intertwined. Whether 
AGNs play a direct role in the evolution of galaxies, 
for example, regulating star formation via feedback pro- 
cesses, is a vital, unresolved question. 

In this section we investigate the positions of AGNs 



in the optical color-magnitude diagram and determine 
the prevalence and distribution of AGN activity within 
the two prominent galaxy populations — blue cloud galax- 
ies, associated with ongoing star formation activity; and 
the red sequence of generally older, quiescent, and more 
massive galaxies — as well as the green valley that lies 
in b etween and may represent a transition population 
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(e-g- JMartin et al.ll2007HSalim et all2009t iMendez et all 
120111 ). Indeed, an overdensity of AGNs in the luminous 
galaxie s close to and within the green valley has been ob- 
serve d (jNandraet al.ll2007t ICoil et ah 2009; H ickox et all 
2009) and may be evidence that AGN activity is associ- 
ated wit h the color trans formation of galaxies. Recent 
work by iXue et al.l (|2010f) . however, has shown the im- 
portance of stellar-mass selection effects in such analyses, 
indicating that stellar mass, not color, may be the key 
parameter that drives the observed trends in the color- 
magnitude diagram. In Sections [5] and [6] above we show 
that p{Lx | A4*, z) is strongly dependent on stellar mass, 
redshift, and X-ray luminosity (although the underlying 
trend is likely driven by Eddington ratio, see Section [7]). 
These trends will significantly affect the observed distri- 
bution of AGNs in the color-magnitude diagram. Thus, 
we must account for these effects to investigate the true 
connection between the prevalence of AGN activity and 
the colors of potential host galaxies, allowing us to shed 
light on the relationship between SMBH accretion and 
star formation processes. 

8.1. AGNs in the Color-magnitude Diagram 

Figure O (top) shows the color-magnitude diagram (af- 
ter applying the redshift-dependent stellar-mass cut de- 
scribed in Section [472] and shown in Figure[3]) for our par- 
ent sample of PRIMUS galaxies and our X-ray-detected 
AGN samples. The color distributions of all galaxies 
change over the redshift range spanned by our sample; at 
higher redshift the position of t he red sequence a nd blue 
cloud shift to bluer colors (e.g.. iBell etai] 120041 ). Addi- 
tionally, the green valley that divides these populations 
has a weak magnitude dependence, following the slope of 
the red sequence (as can be seen in Figure [21 top). Stel- 
lar mass limits also cut diagonally t hrough the color- 
magnitude diagram (|Bell et al.ll2003f) . as demonstrated 
by the blue dashed line in Figure M (top) . To cleanly 
separate blue and red galaxies and reveal any possible 
overdensity of AGNs within the intermediate population, 
we choose to present results based on the color relative 
to the green valley, C, where we define the green valley 
as a function of M g and redshift using EquationQ] Thus, 

C={u-g)- (0.671 - 0.031M g - 0.065z) (18) 

where (u — g) is the rest-frame color of a galaxy and M g 
and z arc the absolute g-band magnitude and redshift. 
We re-project the color-magnitude diagram to the more 
informative C versus A4*, as shown in the bottom panel 
of Figure [9J It is immediately apparent that our detected 
AGN sample is predominantly in the more massive galax- 
ies, consistent with the increase in p(Lx | M*,z) found 
in the preceding sections. They do, however, appear to 
span the full range of galaxy colors at any given stellar 
mass. 

8.2. The fraction of AGN as a function of host galaxy 

color 
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Fig. 9. — Top: color— magnitude diagram for our parent sam- 
ple of galaxies (small gray dots and contours) and X-ray-detected 
AGN sample (colors and symbols correspond to the mass and lu- 
minosity bins used in Figures [4] and |6j . The green valley evolves 
over the redshift spanned by our sample: the dashed line shows 
the center of the green valley at the center of our redshift range, 
z = 0.6 (from Equation [TJ ; the dotted lines show the green valley 
at z = 0.2 and z = 1.0. The light blue dot-dashed line shows the 
track in color-magnitude space based on the simple linear relations 
betwe en mass-to-light ratios and (u — g) color given by Bell ct al. 
(2003) for log M*/A4q = 10.5. Our mass estimates, which use 
multiple photometric bands and redshift information, do not di- 
rectly follow this linear relation, and thus cut through the color- 
magnitude diagram in a complex fashion. Bottom: as in the top 
panel, but converted to the rest-frame (u — g) color relative to the 
green valley (for the magnitude and redshift of a given source) on 
the y-axis and A4* on the ai-axis. 



In the top panels of Figure [TU] we show histograms 
of C for our parent galaxy samples (dashed gray lines) 
in five stellar mass bins, compared to the histograms for 
the X-ray-detected AGN sub-sample (solid colored lines). 
In two of our stellar mass bins we find that the color 
distributions of the AGN host galaxies and the overall 
galaxy population are significantly different based on a 
Kolmogorov-Smirnov (KS) test, which gives probabili- 
ties of 0.79% and 0.00% that the AGN and galaxy sam- 
ples are drawn from the same color distributions for the 
10.5 < logM*/Mo < 11.0 and 11.0 < logM*/M Q < 
11.5 stellar mass bins, respectively. The colors of the 
AGN host galaxies are in general bluer than the overall 
galaxy population. In the other stellar mass bins the K-S 
test finds that there is no significant evidence (at the 99% 
confidence level) for any difference in the color distribu- 
tions, although this may be partly due to the smaller 
samples of AGNs. 

The middle panels of Figure [TU1 show the observed frac- 
tion of galaxies with an X-ray-detected AGN in three 
color bins corresponding to the blue cloud, green valley 
and red sequence. The strong increase in the observed 
fraction of AGNs as A4* increases dominates the behav- 
ior resulting in an overall increase in the AGN fraction 
for all colors at higher stellar masses, but we also see evi- 
dence for an increased AGN fraction in bluer host galax- 
ies. 

The bottom panels show the observed fraction of galax- 
ies hosting an AGN relative to the fraction predicted by 
our model (with no color dependence) for p(Lx | M.*, z) 
determined in Section [51 This accounts for the strong 
increase in p(Lx | -M*, z) with stellar mass and redshift 
and corrects for the effects of the X-ray sensitivity, al- 
lowing us to compare the model with the observed data 
in different color bins and identify deviations that de- 
pend on color. In our lowest stellar mass bin we find 
no evidence for a deviation from the best-fit relation, 
but in the remaining bins there is evidence for a weak 
enhancement of AGN activity in the blue cloud and 
green valley compared to the red sequence, by a fac- 
tor ~ 2 — 4. There are hints that the enhancement is 
greater in the green valley compared to the blue cloud 
for 10.0 < log M*/ 'M & < 11.0 (blue and green points), 
whereas at 11.0 < \ogM.*/M® < 11.5 (yellow points) 
the enhancement is greater in the blue cloud. Given 
the uncertainties in individual measurements, however, 
there is not significant evidence of a difference between 
the enhancement factor for blue cloud and green val- 
ley galaxies, or for a joint 7W* and color dependence. 
In our lowest and highest stellar mass bins there are 
too few objects to identify any color dependence. We 
also note that due to our redshift-dependent stellar mass 
limits (see Wection 14. 2p our lower stellar mass samples 
(log M.*l Mq < 10.5) are restricted to the lower end of 
our redshift range (z < 0.6). Redshift evolution of the 
overall AGN fraction (as determined in Section [S]) is ac- 
counted for in the lower panels of Figure [101 but the 
results could be biased if the relative enhancement in 
the blue cloud and green valley evolves differently to the 
total population. 

8.3. Maximum-likelihood Fitting in Multiple Color Bins 

To quantify these trends, we have repeated our 
maximum-likelihood fitting (Section l6.1[) to constrain the 
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Fig. 10. — Top panels: histograms of the color (relative to the green valley) of the X-ray-detected AGN samples (colored lines) and 
parent galaxy sample (gray dashed lines) in our five stellar mass bins. The histograms are normalized to show the fraction of sources within 
a given population in each color bin. In each panel we give the total number of galaxies and X-ray-selected AGNs in that stellar mass bin; 
we also show the probability, based on the Kolmogorov-Smirnov test, that the two samples are drawn from the same distribution. In the 
highest two Ait bins we find significant evidence that the underlying distributions are different — AGNs appear to be associated with bluer 
host galaxies. For the lower A4* bins we cannot exclude the possibility that the AGNs are drawn from the same distribution of colors 
as the parent galaxy sample. Middle panels: the observed fraction of parent galaxies found to host X-ray AGNs (i.e., the ratio of source 
numbers used in the histograms in the top panels) in three color bins corresponding to the blue cloud, green valley, and red sequence. No 
correction is performed to account for the X-ray incompleteness, redshift evolution, or dependence of the AGN fraction on Lx or Ait- 
Bottom panels: the ratio of the observed fraction of galaxies hosting AGNs to that expected based on our model for p(Lx | A4*,z) and 
accounting for our X-ray incompleteness. Wc find that the observed fraction of AGNs is higher than expected from our model in the green 
valley and blue cloud, and lower than expected in the red sequence. 

shift can be well described by the functional form derived 
in Section [SJ The bottom panel of Figure QTJ however, 
shows that our best-fit normalizations are significantly 
higher for the blue cloud and green valley than for the 
red sequence. This is consistent with the behavior found 
in Section 18.21 above and seen in Figure [TJJJ The stel- 
lar mass dependence drives the trends seen in the color- 
magnitude diagram, and thus more AGN are found in 
massive galaxies, which tend to have redder colors, but 
once this effect is accounted for we have shown that AGN 
activity is enhanced in galaxies with bluer colors associ- 
ated with increased star formation activity. 

8.4. AGN contamination of the host galaxy light 

We note that the optical light of some of our X-ray 
AGN host galaxies may be contaminated by the AGN, 
which could make the observed SED appear bluer than 
that of the actual galaxy. In these cases our derived 
stellar masses should still be robust as the erroneous 
addition of a young stellar component to account for 



functional form of p(Lx | M*, z, C), which describes the 
probability an AGN of luminosity Lx being found in a 
galaxy with a given stellar mass, redshift and color. We 
assume the same power-law dependence on Lx-, -M.* and 
z given by Equation [T"2l Figure HT1 shows the best- fit pa- 
rameters in the three bins of galaxy color corresponding 
to the blue cloud, green valley and red sequence, com- 
pared to the overall best fit from Section 16.31 There 
is weak evidence that the three power-law slopes (7.M) 
7l and j z ) are dependent on galaxy color. The slope 
controlling the stellar mass dependence (tx) may be 
steeper in the blue cloud, and the luminosity dependence 
may have a flatter slope (71,), indicating a higher preva- 
lence of AGNs in blue galaxies with high stellar masses 
and a tendency for blue galaxies to host more luminous 
AGNs. Nevertheless, all the slopes are consistent with 
our overall best-fit relation to within < 3ct. Thus, while 
there are hints of a more complex behavior, the shape 
of the overall trends in the prevalence of AGN activity 
as a function of stellar mass, X-ray luminosity and red- 
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Fig. 11. — Best-fit parameters from the maximum-likelihood fit- 
ting for p(-Lx I Al* , z) in three color bins corresponding to the blue 
cloud, green valley, and red sequence, compared to overall best fit 
parameters and their uncertainties from Section|6](horizontal black 
dashed lines and gray regions). There is weak evidence that the 
parameters describing the power-law dependences on Lx, Af*and 
z depend on the galaxy color. We find very strong evidence that 
the overall normalization is higher for the blue cloud and green val- 
ley than red sequence by a factor of ~ 2. The prevalence of AGN 
activity appears to be mildly enhanced in galaxies with increased 
star formation activity. 

the blue light will contribute relatively little additional 
stellar mass. However, AGN contamination could skew 
the observed color distribution and contribute to the ob- 
served enhancement of AGN activity in the blue cloud 
and green valley. In the most extreme case, an unob- 
scured AGN that contributes a fraction > 30% of the B- 
band optical light in a red, elliptical host galaxy would 
skew the observed colors by up to A(u — g) ~ 1.0, eas- 



ily moving the galaxy fro m the red sequenc e into the 
green valley or blue cloud (|Pierce et a l. 2010a). We have 
identified and excluded sources that exhibit broad emis- 
sion lines in their optical (PRIMUS) spectra, which will 
exclude the most severely contaminated objects where 
an unobscured AGN is clearly dominating over the host 
galaxy optical light. We also set an upper limit on the X- 
ray luminosity of Lx = 10 44 erg s~ 4 , which prevents the 
inclusion of luminous AGNs that are also likely to have 
biased observed colors. Nevertheless, AGN contamina- 
tion may still a ffect the obs e rved co lors of our sample. 
Recent work by IPierce et al.l (|201 0b) comparing the nu- 
clear and outer colors of X-ray-selected AGN hosts at 
z ~ 1 showed that in general the integrated colors were 
unaffected, except in rare cases (< 10% of galaxies) that 
are X-ray bright (Lx > 2 x 10 43 erg s _1 ) and appear un- 
absorbed at X-ray wavelengths. Thus, AGN contamina- 
tion is not expected to severely bias our overall results, 
but could affect our highest luminosity bins. To check 
the extent of any bias we artificially redden the colors of 
all X-ray-detected AGNs in our highest luminosity bin 
(L x > 10 43 - 5 erg s" 1 ) and in the blue cloud (C < -0.1) 
by A(u — g) = 1.0, and repeat our maximum-likelihood 
fitting. We find that the jl, jm, and 7z slopes are al- 
tered but remain within their statistical error bars. The 
overall factor ~ 2 enhancement in the green valley and 
blue cloud is not affected. Thus, any systematic bias will 
not affect our conclusions even in this extreme case of 
severe AGN contamination for all of our most luminous 



8.5. Eddington ratio distribution for different host 
galaxy colors 

In Figure [T2l we show the distribution of Eddington ra- 
tios for our three color bins, again using the N b a /N m d\ 
method to compare the observed fraction with that pre- 
dicted by our best-fit relation determined in Section [7J 
The results also show an enhancement of p(AEdd | M.*,z) 
in the green valley and blue cloud compared to the red se- 
quence (significant at > 3c when combining over the full 
-^Edd range). However, the shape of the distribution of 
Eddington ratios is roughly the same for all colors, with 
approximately the same slope and AGNs found through- 
out the full observable range of AEdd ~ 10~ 4 — 10 _1 . 
AGN activity is clearly enhanced in galaxies that tend 
to have bluer colors (indicating active star formation is 
taking place), but the uniform distribution of Eddington 
ratios is strong evidence that the basic physical processes 
that trigger and fuel the AGN are essentially the same. It 
should be noted that the color-dependent enhancement 
is relatively weak (factor ~ 2) compared to the overall 
trend with Eddington ratio or the strong redshift depen- 
dence. We discuss our findings and their implications 
further in Sections 19.41 and 19.51 below. 



9. DISCUSSION 

The results presented in this paper reveal how the 
prevalence and distribution of moderately obscured, 
moderate-luminosity AGN accretion activity is related 
to the properties of their potential host galaxies during a 
key epoch in the lifetime of the universe, 0.2 < z < 1.0. 
We find that the distribution of X-ray luminosities for all 
stellar masses has the same power-law shape. The prob- 
ability of a galaxy hosting an AGN above a given X-ray 
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Fig. 12. — Distribution of Eddington ratios, p(Asdd I M*, z, C), 
for galaxies in the blue cloud, green valley, and red sequence (top). 
The bottom panel shows the ratio of the data in the top panel to our 
best-fit model from Section [7] We find a slightly higher fraction of 
galaxies hosting AGNs in the blue cloud and green valley (factor ~ 
2), although the distribution of Eddington ratios has approximately 
the same shape, with AGNs found over the full range of Eddington 
ratios in galaxies of all colors. 

luminosity increases strongly with the stellar mass of the 
galaxy. However, we show that this effect may be pri- 
marily driven by the Eddington ratio distribution. The 
distribution of Eddington ratios has a power-law shape; 
a low fraction of galaxies (~ 0.1%) host AGNs accreting 
close to the Eddington limit (AEdd<; 0.1), whereas ~ 10% 
of galaxies contain an SMBH growing at a very low rate 
(AEdd~ 10~ 4 ). The probability of hosting an AGN of a 
given Eddington ratio is independent of stellar mass. 

Alternatively, we can interpret these results in terms 
of a specific accretion rate: the rate of black hole growth 
relative to the stellar mass of the galaxy. The distri- 
bution of specific accretion rates is a universal function 
that is independent of stellar mass (see Section I7H|) . Our 
conclusions hold regardless of uncertainties in the black 
hole mass and thus the absolute value of the Eddington 
ratios for the AGN. 

We find that AGNs are not predominantly in massive, 
red galaxies. More AGNs may be found in red galaxies 
due to a combination of two selection effects: (1) it is 
easier to detect the more prevalent, low Eddington ratio 
(low specific accretion rate) AGNs in galaxies with higher 
stellar masses as the AGNs will be intrinsically brighter; 
and (2) more massive galaxies tend to have redder colors. 

There is, however, a mild dependence of the preva- 
lence of AGN activity on galaxy color. While AGNs are 
found in galaxies across the whole range of colors, there 
appears to be an enhancement (factor ~ 2) in the prob- 
ability of hosting an AGN for galaxies with blue colors, 



both within the blue cloud and the green valley. Thus 
there may be a weak association between AGN activity 
and star formation within a host galaxy. Nevertheless, a 
lack of star formation does not preclude AGN activity — 
the distribution of Eddington ratios in red galaxies fol- 
lows approximately the same power-law shape as in blue 
galaxies, spanning the full range of accretion rates. 

AGN activity also evolves strongly with redshift. We 
find that the probability of a galaxy hosting an AGN has 
decreased rapidly since z ~ 1 , but the distribution of Ed- 
dington ratios retains the same power-law shape. Thus, 
the same physical process or processes are likely to be 
responsible for triggering and fueling AGNs throughout 
this epoch, but they must diminish rapidly over the time 
period (see Section [973]) . 

In this section we discuss the implications of our re- 
sults. In Sections 19.11 and 19.21 we compare our findings 
with previous studies of the Eddington ratio distribution 
at similar redshifts to our study and in the local universe 
respectively. In Section [9731 we discuss the implications 
of our results for models of AGN triggering and fuel- 
ing processes over the last 8 billion years. In Section l9~4l 
we discuss our findings on the dependence of AGN accre- 
tion activity on host-galaxy color and compare with prior 
studies of color-magnitude relations for AGNs. Section 
9.51 discu sses the role of AGNs in galaxy evolution. Sec- 
tion [9?6] discusses the remaining issues and uncertainties 
in our work and prospects for future studies. 

9.1. Comparison with prior studies of the Eddington 
ratio distribution at z > 0.2 

Several prior studies have presented distributions of 
Eddington rati os for X-ray- s electe d samples of AGNs 
at z > 0.2. iHickox et all (|2009l ) found that X-ray- 
selected AGNs in AGES at 0.25 < z < 0.8 have a wide 
range of Eddington ratios spanning 10~ 3 < AEdd < 1, 
with a broad, roughly lognormal distribution peaking 
at AEdd ~ 0.02. However, the incompleteness effects 
due to the X-ray sensitivity were not taken into ac- 
count. iSimmons et all (2011) determined Eddington ra- 
tios for X-ray sources in the GOODS fields (using very 
deep X-ray data from both the Chandra Deep Fields) at 
0.25 < z < 1.25 and found that the sample was dom- 
inated by objects accreting at very low Eddington ra- 
tios (AEdd < 0.01) with a median of Asdd = 0.006; the 
extremely deep X-ray data mitigate incompleteness ef- 
fects an d enable detection o f very low Eddington ratio 
sources. IBrusa et all (|2009bD also presented an Edding- 
ton ratio distribution for obscured, moderate-luminosity, 
X-ray AGNs identified in the Chandra Deep Field-South 
at higher redshifts: 1 < z < 2 and 2 < z < 4. Simi- 
lar, broad distributions of Eddington ratios were found 
in both redshift bins with medians AEdd ~ 0.02 — 0.08. 
None of these prior papers have tracked or corrected for 
X-ray selection effects to reveal the underlying distribu- 
tion of Edd ington ratios. 

Recently, iTrump et all (|2011l ) proposed that there are 
two distinct distributions of Eddington ratios for X-ray- 
selected AGNs corresponding to two different modes of 
accretion. They find that broad-line QSOs are predomi- 
nantly at higher AEdd~ 0.1, consistent with other stud- 
ies of the Eddington ratio distribution of QSOs (e.g. 



Kollmeier et"alll2006l: iGavignaud et~alll2008l : iKellv et all 
2010T) . Thev propose that these AGNs are fed by a thin 
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accretion disk with a disk wind containing the broad-line 
region and some obscuring material. Narrowdine AGNs 
or X-ray AGNs lacking optical AGN signatures are shown 
to have lower AEdd~ 10 -3 and may be powered by a ge- 
ometrically thick, radiatively inefficient accretion flow. 
Our sample was explicitly limited to "obscured" AGNs 
that lack broad lines in their optical spectra (allowing us 
to probe the host galaxy properties in the optical light) 
and will thus corres pond to the second class proposed by 
I Trump et al.l (J2011I ). Consistent with their findings, we 
have shown that these AGNs are predominantly at low 
Eddington ratios. In fact, we find a power-law relation 
that extends to very low AEdd^ 10 -4 . Therefore, the 
physical processes that trigger and fuel low-level AGN 
accretion activity may also be common and take place in 
galaxies of all stellar masses. 

Our results are consistent with previous studies of 
Eddington ratio distributions for moderately obscured, 
moderate-luminosity AGNs at z > 0.2. However, we 
show that the true distribution of Eddington ratios is a 
power law that continues to rise to lower values of AEdd- 
This result is reliant on careful analysis of X-ray sen- 
sitivity effects relative to the stellar masses of the host 
galaxies. Our directly observed distribution of Edding- 
ton ratios has a broad, approximately lognormal distri- 
bution with a peak at AEdd^ 10~ 2 (see gray histogram in 
final panel of Figure [7]) as found in prior works but this is 
significantly skewed by X-ray incompleteness effects and 
does not represent the true AEdd distribution. We note 
that our study has excluded BLAGNs which may domi- 
nate the population at high AEdd and could correspond 
to a distinct accretion mode. 

9.2. Comparison with prior studies of the local 
Eddington ratio distribution 

The lack of an extremely large area, sensitive X- 
ray survey at > 2 keV has meant that the major- 
ity of studies of the Eddington ratios of AGNs in 
the local universe {z < 0.1) have relied on large-area 
optical spectroscopic surveys to identify AGNs from 
their optical line-emission, in particular the Sloan Dig- 
ital Sky Survey (SPSS: e.g.. IKauffmann et al.l 12 003a: 
Heckman et al.l'l2004j: IKauffmann et al.l 12004 lHao et al.l 
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though see also iGeorgakakis fc Nandrall20lTj) . As dis- 
cussed in Section [1] such studies have found that the 
fraction of galaxies hosting optically selected narrow- 
emission-line (Type 2) AGNs rises steeply with increas- 
ing stellar mass and that Type 2 AGNs are found 
predominantly in massive galaxies (M. * > 10 10 Mq\ 
IKauffmann et al.ll2003aHBest et al.l l2005ft. 

Heckman ct al. (2004, hereafter H04) studied the 
present-day growth rates of black holes in SDSS galaxies 
using the Oin luminosity as a tracer of AGN accretion ac- 
tivity. They found that the distribution of Eddington ra- 
tios is highly dependent on black hole mass; while ~ 0.5% 
of black holes with .Mbh = 10 7 A / fo are accreting at Ed- 
dington ratios AEdd > 0.1, less than 0.01% of higher mass 
black holes (A^bh ~3x 10 8 .M^) are accreting a t these 
high Eddington ratios. IKauffmann &; Heckman! (2009, 
hereafter K09) built on this work and proposed that black 
hole growth in the local universe is driven by two distinct 
processes: (1) a self-regulated mode of growth in young 
galaxies with plentiful cold gas that leads to a broad 



lognormal distribution of Eddington ratios peaking at 
AEdd^ 10 -2 - 5 , and (2) a passive mode in galaxies with 
older central stellar populations where the supply of fuel 
to the black hole is driven by stellar mass loss from stars, 
resulting in a power-law distribution of Eddington ratios. 
This is somewhat different to our results at higher red- 
shifts using X-ray selection but probing similar ranges of 
Eddington ratios and black hole masses (roughly traced 
by the galaxy stellar mass). We find that the distribu- 
tion of Eddington ratios is the same for all stellar masses 
and is described by a power law; we find no evidence for 
an additional lognormal component and that the shape 
of the Eddington ratio distribution is similar for young 
and old stellar populations (i.e., blue and red galaxies). 

It is unclear how to reconcile our findings with the re- 
sults of H04 and K09 without further study to directly 
compare X-ray and optical probes of AGN accretion in 
both the local and high-redshift universe. The discrep- 
ancies could be related to biases inherent to the two very 
different selection techniques that affect both the com- 
pleteness of the AGN samples and the estimates of the 
Eddington ratios. For example, AGNs in the lognormal 
mode identified by K09 could be heavily obscured and 
thus are not detected in X-rays or are detected but their 
luminosities (and thus Eddington ratios) are significantly 
underestimated. Alternatively, it is possible that the Oni 
emission line may not provide an accurate tracer of the 
black hole accretion at low Eddington ratios or in galax- 
ies with current star formation. 

Another possibility is that the differences between our 
results and those of H04 and K09 are due to evolution 
with redshift. Indeed, H04 consider volume-averaged 
growth rates and show that the most massive black holes 
in their sample are not growing at significant rates at the 
present day and must have undergone periods of rapid 
growth in the earlier universe. At higher redshifts we may 
be seeing this earlier phase of relatively rapid growth for 
all galaxies over the range of stellar masses (and thus 
black hole masses) that we are able to probe, although 
our results do not show any direct evidence for such 
"downsizing" behavior (see further discussion in Section 
EH below). 

The relative importance of the two accretion modes 
proposed by K09 could also change with redshift. We find 
that the distribution of Eddington ratios has a power-law 
shape with a normalization that increases rapidly with 
redshift as ~ (1 + z) 35 . The lognormal accretion mode 
seen by K09 may not evolve or may evolve only weakly 
with redshift and thus would not be seen in our sample. 
This may, however, conflict with the K09 interpretation 
that the accretion is driven by stochastic accretion of cold 
gas, as the supply of cold gas should be more plentiful in 
the high-redshift universe. 

Further work and a direct comparison of X-ray and 
optical probes of AGN accretion modes in both the local 
and high-redshift universe is required to resolve these 
issues. 

9.3. The Evolution of AGN Triggering and Fueling 

Our observed AGN Eddington ratio distribution can 
provide key insights into the physical processes that 
are responsible for triggering and fueling AGNs at z ~ 
0.2 — 1. One of our major findings is that the probability 
of hosting an AGN for galaxies throughout our stellar 
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mass range (9.5 < log M^/Mq < 12) is primarily deter- 
mined by the Eddington ratio distribution and is inde- 
pendent of the stellar mass. Not only is the shape of the 
distribution of accretion rates (tracked by the Eddington 
ratios) the same, the overall normalization does not de- 
pend on the stellar mass. Thus, AGN activity may take 
place in any galaxy, regardless of the stellar mass. The 
universality in both the prevalence of black hole growth 
and the distribution of accretion rates (Eddington ratios) 
indicates that the same underlying physical processes are 
responsible for triggering and fueling AGN activity in all 
moderately massive galaxies. 

We also find that the fraction of galaxies hosting AGNs 
above a given Eddington ratio evolves strongly with red- 
shift with the form (1 + 2;) 3 - 5±0 - 5 to z = 1. The shape 
of the distribution remains the same, indicating that the 
same processes are regulating AGN activity but are un- 
dergoing a rapid evolution with redshift. This evolution 
could take two possible forms: (1) the fraction of ac- 
creting systems is higher at higher redshift, or (2) the 
characteristic Eddington ratio is higher at higher red- 
shift. In both scenarios, a given physical process must 
trigger AGN activity across our redshift range, resulting 
in the same distribution of Eddington ratios with time. 
However, in the first scenario the fraction of galaxies with 
this process is higher at earlier times, while in the second 
scenario the AGNs are more rapidly accreting at earlier 
times. We cannot distinguish between these two scenar- 
ios solely with the results presented here. 

However, studies of the evolution of the X-ray lumi- 
nosity function (XLF) of AGNs, which traces the overall 
distribution of AGN accretion activity, find that at low 
redshifts {z < 1) the evolution is dominated by lum inos- 
ity evolution (jBarger et alJl2005t lAird et alll2010D . The 
characteristic luminosity, L» , corresponding to the break 
in the double-power-law form of the XLF, shifts to lower 
luminosities between z ~ 1 and the present day. The 
Eddington ratio distribution presented here will track 
the evolution of the faint end of the XLF. Indeed, if 
we translate the luminosity evolution of L» into an ef- 
fective change in normalization at the faint end of the 
XLF w e find evolution o f the form (1 + z ) 39±0 - 3 (based 
on the lAird et aLl|2010i "luminosity and density evolu- 
tion" model), in excellent agreement with our evolution 
in the normalization of the Eddington ratio distribution 
presented above. We therefore conclude that our redshift 
evolution is the result of a shift in the accretion processes 
to lower characteristic Eddington ratios since z ~ 1. 

This implies that any upper boundary to the Edding- 
ton ratio distribution would also move to lower AEdd as 
redshift decreases. Our results do not show evidence for 
such a redshift-dependent cutoff, although we probe only 
relatively low Eddington ratios, especially at low red- 
shifts (see Figured]), and do not include broad-line QSOs 
that may dominate the population at the highest Edding- 
ton ratios. Much larger samples are required to deter- 
mine the presence and evolution of an upper boundary 
or turnover o f our E ddingt on rat i o dist ribution at high 
AEdd- In fact. lSteinhardt fe Elvisl (|2010l ) have presented 
evidence for a "sub-Eddington boundary" in large sam- 
ples of SDSS quasars at 0.2 < z < 4, which depends on 
black hole mass and may evolve with redshift, although 
the connection to our results requires further investiga- 
tion. 



Our findings do not agree with the "downsizing" pic- 
ture of AGN evolution, which attributes the luminos- 
ity evolution of the XLF to a decrease in the aver- 
age mass of galaxies that host AGNs as time pro- 
gresses, either due to a shut off of AGN activity in 
the highest mass galaxies or AGNs i n lower mass galax- 
ies "turning on" at later times (e.g.. iBarger et alj|2005t 
lAlonso-Herrero et all 120081; IHopkins fe Hernquistl 12009). 
This picture is not supported by our results: we find that 
the distribution of accretion rates is uniform for galaxies 
of all stellar masses at any given redshift. The entire dis- 
tribution evolves with redshift, which we interpret as a 
shift to lower characteristic accretion rates with increas- 
ing cosmic time. Thus, at any given epoch, black hole 
growth is taking place in galaxies over the full range of 
stellar masses, but at later times the accretion rates are 
simply much lower, on average, and there is a relative 
lack of the most rapidly accreting sources. The evolu- 
tion of the XLF is driven by this shift in accretion rate 
rather than the dow nsizing of black hole growth (see also 
Shan kar et all l200"9f ) . However, as our results track the 
specific accretion rate relative to the host stellar mass 
(see Section l7.4|) . we are unable to constrain any under- 
lying evolution of the distribution of black hole masses, 
which may evolve differently from the stellar masses of 
their host galaxies and could potentially exhibit a down- 
sizing behavior. 

While our results provide important constraints, it is 
currently unclear what physical mechanism could be re- 
sponsible for fueling AGN activity and regulating the 
distribution and evolution of accretion rates. The ob- 
served drop in the global star formation rate density 
from z ~ 1 - 2 to z ~ (Hopkins I2004D presumably 
reflects a drop in the availability of cold gas in galax- 
ies to fuel star formation and potentially AGNs. How- 
ever, gas must be driven to the central region of the host 
galaxy to fuel an AGN. The fueling of the AGN could be 
triggered by instabilities i n the host galaxy caused eithe r 
by secular evolution (e.g.. lKormendv fc Kennicuttll2004|) . 
galaxy-gal axy interactions, or major o r minor merg- 
ers ( e.g., Kauffma nn fc Haehneltl 120001 : IHopkins et all 
12008( 1. Alternatively AGN activity could be fed 
by smooth cosm ological cold gas accretion (e.g., 
iKeres etlrtl 120051). stochastic accretion of g as within 
the galaxy (e.g., IHopkins fc Hernquistl 120061) . or mass 
loss from evolved stars (e.g ., ICiotti fc Ostrikerl 120071 : 
iKauffmann fc H cckman 2009). Whichever process con- 
trols the fueling of AGNs must be relatively independent 
of the stellar mass of the host galaxy, at least over the 
range probed by our study. Significant additional study 
is required to fully reveal the processes that are driving 
AGN accretion activity. 

9.4. The relationship between AGN activity and host 
galaxy color 

We find that the prevalence of AGN activity is some- 
what enhanced in galaxies with blue optical colors, 
once stellar-mass-dependent selection effects are fully ac- 
counted for. In addition, we show that the probability 
of hosting an AGN for galaxies throughout our stellar 
mass range is primarily determined by the Eddington 
ratio. This drives the observed increase in the frac- 
tion of galaxies hosting X-ray- detected AGNs as the stel- 
lar masses of the parent galaxy sample increase. Thus, 
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galaxies throughout the color-magnitude diagram may 
host an AGN, but a higher fraction of AGNs are detected 
within the most luminous galaxies, which also tend to 
have redder colors. In fact, the probability of hosting an 
AGN above a given Eddington ratio is higher for galaxies 
within the blue cloud and the green valley. 

Both Kue et al.l (pOloT) and (Silverman et all (|2Q09bf ) 
have previously investigated the fraction of AGNs as a 
function of host - galaxy color and stellar mass. Indeed, 
ISilverman et a l. (2009b) highlighted the importance of 
stellar-mass-dependent selection effects and the potential 
bias in studies using optical flux-limited samples (e.g. , 
Nandra et all 120071 iCoil et al.ll200i IHickox et all [20091: 



Schawinski et al.l l2009t iGeo rgakaki s fc Nandral |20TTT) 
which was la ter demons t rated explicitly by Kue et al] 
(2010). The Kue et al.l (|2010D study also found that 
moderate-luminosity X-ray AGNs are distributed over 
the full range of colors in stellar-mass-matched samples 
with no specific clustering in the color-magnitude dia- 
gram. However, their sample of AGNs at < z < 1 is 
comparatively small (83 X-ray-detected AGNs) and was 
combined over this entire redshift range where we have 
found strong evolution of the AGN fraction. In addi- 
tion, they are mainly reliant on photometric redshift esti- 
mates, which have larger errors and a higher outlier frac- 
tion than our PRIMUS redshifts, and probe to fainter op- 
tical magnitudes. These uncertainties will propagate to 
estimates of stellar masses and rest-frame colors, poten- 
tially blurring the color-magnitude diagram and mask- 
ing the relatively weak enhancement (factor ~ 2) of the 
AGN fraction that we see in galaxies in the blue cloud 
a nd green valley. 

ISilverman et all (|2009bl ) found that the X-ray AGN 
fraction at 0.5 < 2 < 1.0 increases in galaxies with young 
stellar ages, traced by the D„(4000) spectral index, and 
in galaxies with blue rest-frame optical co lors, consistent 
with our results. ISilverman et al.l (|2009bl ) also presented 
evidence that the levels of ongoing star formation, traced 
by the [On]A3727 emission-line luminosity, are higher in 
AGN host galaxies than in non-AGN galaxies with equiv- 
al ent stellar masses . 

ISchawinski et al.l (|2010t herafter S10) also studied the 
AGN fraction as a function of galaxy color and stellar 
mass. They used a local sample of galaxies from the 
SDSS (z < 0.05), identified narrow-line AGNs within 
these galaxies using emission-line diagnostics and deter- 
mined AGN luminosities (and thus Eddington ratios) 
from the intensity of the [OI 1 1] emission line. S10 find 
that the fraction of detected AGNs increases with stel- 
lar mass for all galaxies, consistent with prior studies 
and our findings for X-ray-detected AGNs. An enhance- 
ment in the fraction of AGNs is seen in the green val- 
ley along with a significant lack of emission-line AGNs 
in galaxies on the red sequence. They also find a low 
fraction of AGNs in the blue cloud but note that for blue 
galaxies their AGN sample is incomplete for Loin < 10 41 
erg s _1 ; at low luminosities the AGN emission lines are 
overwhelmed by the star formation precluding identifica- 
tion of the AGN. The blue galaxy sample is also gener- 
ally at low stellar masses (log A'U/.M© < 10.5). Using 
our best-fit results for p(AEdd | M*,z) from Section [7] 
and assuming a bolometr i c corr ection of Lboi ~ 450Loiii 
(|Kauffmann fc Heckmar] I2009D . we predict < 1% of 
galaxies with 9.5 < logM*/M Q < 10.5 at z < 0.05 



to host AGN with L iii > 10 41 . Thus the results of S10 
appear consistent with our findings, but the bias against 
the selection of AGNs in low-mass, star-forming galaxies 
in their work may have a significant effect. 

9.5. The role of AGNs in galaxy evolution 

What do our results imply about the connection be- 
tween AGN activity and galaxy evolution? The frac- 
tion of galaxies hosting AGN rapidly decreases between 
z « 1 and z ss 0.2 but follows the same distribution of 
Eddington ratios. We find strong evolution in the frac- 
tion of galaxies hosting AGN with redshift of the form 
p(AEdd | M*,z) oc (1 + z ) 3 - 5±0 - 5 . This evolution is very 
similar to the overall evolution in the s tar formation den- 
sity, p „ ex (l+z) 35±0 ' 2 for < z < 1.2 (|Ruionakarn et all 
2010). This correspondence could be due to a relation- 
ship between the processes that fuel star formation and 
AGN: both could be fueled from the same gas supply, or 
the AGN could be fueled indirectly via stellar mass loss 
from the stars within the gala xy (e.g. JNorman fc Scovillel 
198*, iCiotti fc Osilikenl20q7l) . 

We also find a weak association between the prevalence 
of AGN activity and blue optical colors associated with 
current star formation activity. That said, AGNs with 
a wide range of Eddington ratios are found in galaxies 
of all colors. The shape of the distribution of Eddington 
ratios is roughly constant as a function of color. The en- 
hancement of AGN activity in galaxies with blue colors 
could indicate that the distribution of Eddington ratios 
is shifted to higher AEdd, on average. This could imply 
that red galaxies contain AGNs at later stages of their 
lifetimes when the rate of accretion has reduced and the 
luminosities have faded. It has been proposed that feed- 
back from an AGN can quench star formation in galax- 
ies, resulting in the re d colors (e.g.,[Springcl et al. 2005; 
iHopkins et al.l l2006bf l . However, the most rapidly ac- 
creting SMBHs (highest Eddington ratio sources) in our 
sample are not exclusively in galaxies with current star 
formation (blue colors). Red colors and a lack of current 
star formation for a galaxy do not preclude substantial 
AGN accretion and significant SMBH growth. Equally, 
a high fraction of blue galaxies contain weakly accreting 
(low Eddington ratio) AGNs. 

Thus, while AGN activity and star formation do ap- 
pear to be correlated, we do not find evidence that AGNs 
play a direct role in shaping the global properties of their 
host galaxies or their evolution. However, galaxies and 
their central AGNs could co-evolve in that the AGN frac- 
tion, related to the duty cycle of AGNs, may be regulated 
by whatever process drives gas to the centers of galax- 
ies and fuels black hole growth. As discussed in Section 
19.31 above, this process must result in a shift to lower 
characteristic accretion rates as cosmic time progresses. 
This does not prevent the triggering of AGN accretion in 
galaxies of any color but must lead to characteristically 
lower accretion rates in red galaxies or a reduction in 
the triggering rate for red galaxies. As with the redshift 
evolution of the Eddington ratio distribution, we cannot 
distinguish between a shift to characteristically higher 
AEdd for blue galaxy colors or an overall increase in the 
density of blue galaxies hosting AGNs at a given AEdd- 

We thus conclude that our results do not show evidence 
that feedback from moderate-luminosity AGNs plays a 
direct role in the evolution of galaxies, although star for- 
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mation and AGN activity may have a common triggering 
mechanism. We note that our results have no bearing on 
any potential role that may be played by luminous QSOs 
in the evolution of galaxies. 

9.6. Remaining issues and future prospects 

Despite our unprecedentedly large sample of galaxies 
and hard X-ray-selected AGNs with spectroscopic red- 
shifts at z ~ 0.2 — 1, our results are dominated by Poisson 
uncertainties. This is due to the relatively small fraction 
of galaxies that host AGNs above a given X-ray lumi- 
nosity (~ 1%— 10% but highly dependent on stellar mass 
and redshift) as well as the depths of our X-ray data 
which impede detection of lower luminosity AGNs. We 
account for the Poissonian uncertainties in our analysis 
and are able to constrain functional forms that provide a 
good description of our data, although our discriminat- 
ing power is limited by the source numbers. Indeed, the 
results of Section[5]hint that the probability of hosting an 
AGN may have a complex color dependence that deviates 
from our global Eddington ratio dependence. However, 
due to our limited sample, we cannot split the data into 
many fine bins in redshift, stellar mass, X-ray luminos- 
ity and galaxy color at the same time to examine this 
behavior more closely. 

We also note that our results are valid for the popula- 
tion of moderately obscured, moderate- luminosity AGNs 
probed by our data. This population appears to rep- 
resent the majority of the AGN accretion density ove r 
our redshift range (|Ueda et al.1 120031 : lAird et aljfeoiOj l. 
and thus our results can shed light on the dominant pro- 
cesses that are responsible for triggering and regulating 
AGN activity over this epoch. Our results should not 
be extrapolated to the very low or high Eddington ratio 
regime where different accretion and fueling processes 
may be relevant. Furthermore, the lowest stellar mass 
galaxies in our sample are only included at the lower end 
of our redshift range due to the redshift-dependent stellar 
mass cut, which is defined by the magnitude limits of our 
survey. This cut could introduce a bias in our results, al- 
though applying an additional cut of log .M*/.M© > 10.5 
at all redshifts does not significantly alter our results (but 
does substantially increase the errors due to the smaller 
sample size). Conversely, higher stellar mass galaxies 
tend to be found at higher redshifts due to the larger vol- 
ume surveyed, which could also introduce a bias. Both 
larger area surveys and deeper samples are required to 
address these issues. 

BLAGNs may dominate the population at the high- 
est Eddington ratios, and thus their exclusion could also 
bias our derived Eddington ratio distribution. To inves- 
tigate the impact of BLAGNs, we apply a completeness 
correction term in our maximum-likelihood fitting based 
on the fraction of BLAGNs within the total sample of 
X-ray detected PRIMUS sources at 0.2 < z < 1.0. We 
apply this correction as a function of X-ray luminosity 
(over the 42 < logix < 44 range) but assume no red- 
shift evolution and implicitly assume that the host stellar 
masses of the BLAGNs have the same distribution as the 
non-broad-line sample. Using this basic correction, we 
find a slightly flatter slope, 7# = 0.62 ± 0.04, which is 
consistent with the slope derived in Section [7] (within the 
errors) and indicates that the exclusion of BLAGNs does 
not introduce a significant bias. Nevertheless, this sim- 



ple correction may not accurately account for this pop- 
ulation. BLAGNs may correspond to a specific phase 
in the lifetimes of AGNs, be powered by different accre- 
tion modes, or be a distinct population that is triggered 
and fueled by very different processes compared to our 
sample. Indeed, BLAGNs dominate the X-ray detected 
population at high luminosities (Lx > 10 44 erg s" 1 ), out- 
side the range probed by our study and corresponding 
to the strong break seen in the XLF. Such sources may 
dominate the AGN population at the highest Eddington 
ratios and constraints on their accretion properties are vi- 
tal to track the shape of the Eddington ratio distribution 
and any cutoff around the Eddington limit. Estimates of 
the SMBH mass for such sources are possible using high- 
resolution spectroscopy and scaling relations between the 
lumi nosity, the width of broad lines and the SMBH mass 
(e.g.. lTrump et afll2011[ ). which would allow estimates of 
the Eddington ratios and a more thorough investigation 
of the Eddington ratio distribution of BLAGNs. How- 
ever, this is not possible with our PRIMUS spectra and 
is thus beyond the scope of this paper. 

Our hard X-ray selection may also bias our sample 
against the most absorbed AGNs (Ah > 10 23 cm~ 2 ) 
and very heavily obscured, Compton-thick sources will 
be missed completely. While such sources may corre- 
spon d to a substantia l fraction of the AGN popu lation 
(e.g.. iGilli et~aL|[2007t iDraper fc Ballantvnd 12009). it is 
unclear if they represent a distinct phase of AGN growth 
and thus how their co-evolution with galaxies may differ 
to the moderately obscured population. 

Our work shows that selection effects can severely bias 
samples of AGNs selected using different methods due to 
the interplay between the Eddington ratio, the observed 
AGN luminosity, and the stellar mass and star formation 
properties of the galaxy. However, by correcting for the 
various selection effects in our data, we are able to reveal 
the overall trends in the prevalence of AGNs and the 
dependence on redshift, stellar mass, X-ray luminosity 
or Eddington ratio. Thus, our results can strongly con- 
strain models of AGN triggering and the co-evolution of 
AGNs and galaxies. Future work could push our studies 
to higher redshifts or compare our findings with rela- 
tionships in the local universe. Alternatively, increased 
area coverage and X-ray depth could increase our sam- 
ple size and allow a more thorough analysis of the rela- 
tionship between luminosity, absorption properties, bolo- 
metric corrections, accretion rates, stellar masses, galaxy 
colors, star formation properties, and morphologies as 
a function of redshift over this key period of the uni- 
verse, providing further insight into the physics of black 
hole growth and the relationship between black holes and 
their host galaxies. However, any such future studies 
must carefully account for stellar-mass-dependent selec- 
tion effects and the X-ray sensitivity to shed light on the 
true distribution of AGN accretion. 

10. SUMMARY 

In this paper we study how the probability of a galaxy 
hosting an AGN depends on stellar mass and color, as 
well as revealing the underlying distribution of accretion 
rates probed by the X-ray luminosities and Eddington ra- 
tios. We use an unprecedentedly large sample of galaxies 
with spectroscopic redshifts from PRIMUS, deep broad- 
band optical imaging, and X-ray data from Chandra or 
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XMM-Newton. 
Our main conclusions are as follows: 

• We find that for galaxies at z = 0.2 — 1 the prob- 
ability of hosting an AGN as a function of X-ray 
luminosity, p(Lx \ M.^ 1 z) 1 is a power-law function 
of M. % with slope "fM that is independent of Lx- 
We also find that the distribution of X-ray lumi- 
nosities at fixed stellar mass is a power law with 
slope 7l that is independent of stellar mass. These 
results imply that while AGN activity is enhanced 
in higher stellar mass sources, the distribution of 
accretion rates is independent of mass. 

• We find that the power-law slopes do not change 
with redshift, but there is strong evolution of the 
overall normalization. The probability of a galaxy 
hosting an AGN drops by a factor ~ 6 between 
z ~ 1 and z ~ 0.2. 

• We show that the probability of a galaxy hosting 
an AGN is defined by a universal Eddington ratio 
distribution, which is a power-law function of AEdd 
with slope 7s = —0.65 and is independent of the 
stellar mass of the host galaxy. 

• Alternatively, these results can be interpreted in 
terms of specific accretion rate — the rate of black 
hole growth relative to the stellar mass of the 
host galaxy — rather than the true Eddington ratio. 
Thus, the distribution of specific accretion rates is 
described by a universal power-law function with 
slope —0.65 for all stellar masses. 

• The Eddington ratio distribution evolves strongly 
with redshift with the form (1 + z ) 3 - 5 ±°- 5 j similar 
to the rapid decrease in the total star formation 
rate density since z ~ 1. However, the drop in the 
overall black hole accretion rate density is driven by 
this decrease in the probability of hosting an AGN 
for galaxies of all stellar masses. We do not find 
a shift in the average stellar mass of a galaxy with 
an accreting black hole to lower masses or a shut 
down of AGN activity in the most massive galaxies 
as cosmic time progresses. 

• We find that the probability of a galaxy hosting 
an AGN is mildly enhanced in galaxies with rel- 
atively blue colors by a factor ~ 2, once stellar- 
mass- and rcdshift-dcpcndcnt selection effects are 
fully accounted for. However, AGN activity is tak- 
ing place in galaxies of all colors. The shape of the 
distribution of AGN accretion rates does not ap- 
pear to significantly vary between the blue cloud, 
green valley and red sequence, and the distributions 
span the full range of Eddington ratios. 

• AGNs are not predominantly in red, massive galax- 
ies. The previously observed increase in the AGN 



fraction with stellar mass is a selection effect that is 
driven by the AEdd distribution — a higher fraction 
of AGNs are detected in high stellar mass galaxies 
(with more massive black holes and generally red- 
der colors) as sources accreting at low Eddington 
rates are intrinsically more luminous than in lower 
stellar mass galaxies. 
Our results imply that the underlying physical pro- 
cesses that are responsible for triggering and fueling ob- 
scured AGN activity do not change between z ~ 1 and 
z ~ 0.2 but must decrease in frequency or shift to charac- 
teristically lower Eddington ratios. The same processes 
must take place in galaxies of all stellar masses and 
all colors, although activity appears to be enhanced in 
galaxies with current star formation. While AGN activ- 
ity and star formation appear to be globally correlated, 
we do not find evidence that AGN feedback plays a di- 
rect and important role in the quenching and regulation 
of star formation in galaxies. 
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